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Abstract 
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I. INTRODUCTION 


Most of states of a nucleus are embedded in the continuum of decay channels due to which 
they get a hnite lifetime. That means: the discrete states of a nucleus shade off into resonance 
states with complex energies {Ei — IT,). Ei give the positions in energy of the resonance 
states while the widths Tj are characteristic of their lifetimes. The E^ may be different 
from the energies of the discrete states, and the widths Tj may be large corresponding to 
a short lifetime. Nevertheless, there is a well dehned relation between the discrete states 
characterizing the closed system, and the resonance states appearing in the open system. 
The main difference in the theoretical description of quantum systems without and with 
coupling to an environment of the decay channels is that the function space of the system is 
supposed to be complete in the hrst case while this is not so in the second case. Accordingly, 
the Hamilton operator is Hermitian in the first case, and the eigenvalues are discrete. The 
resonance states, however, characterize a subsystem described by a non-Hermitian Hamilton 
operator with complex eigenvalues. The function space containing ’everything’ consists, in 
the second case, of system plus environment. 


The mathematical formulation of this problem goes back to Feshbach Q who introduced 
the two subspaces of the Hilbert space containing the (i) discrete and (ii) real-energy scat¬ 
tering states, respectively. This has given foundation of the Continuum Shell Model (CSM) 
j2, |3|,|^. Number of particles in the scattering states provides in CSM a natural hierarchy of 
approximations which, phenomenologically, is associated with a successive opening of more 
and more complicated decay channels at higher excitation energies. The technical difficulties 
in the practical implementation of this strategy are such that all past numerical applications 


ving one-nucleon decay 
y proposed numerical 


of the CSM have been restricted to the description of processes invo 
channels only. Few attempts to treat the multiparticle continuum 
schemes which have never been adopted in realistic multiconhguration mixing calculations. 

Feshbach succeeded in formulating a unihed description of nuclear reactions with both di¬ 
rect processes (the short-time scale) and compound nucleus processes (the long-time scale). 
A unihed description of nuclear structure and nuclear reaction aspects is much more compli¬ 


cated and became possible only at the turn of the last century 


m 


ihe framework of the 


so-called Shell Model Embedded in the Continuum (SMEC) (see Ref. ^ for recent review) 
which is parented to the CSM. The SMEC has been applied for the description of spectra 
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tering continuum, like the (p,p') reaction, 


10|, the Coulomb dissociation reaction 


and reactions involving one particle in the sea 
the proton/neutron radiative capture reactions 

^1, or the hrst-forbidden /3-decay 121. Further applications of the SMEC with one-particle 
continuum involved the study of the continuum effects in the binding systematics of neutron- 


Q. 


and the statistical aspects of the continuum coupling for unbound 


rich nuclei in sd shell 
states in ^^Mg Q. 

First multiconhgurational shell model (SM) approach with no restriction on the number 

in the complete Berggren 


of particles in the continuum has been proposed recently 15, Q ii 

single-particle (s.p.) basis |l^ which consists of Gamow (or resonant) states and the complex 
non-resonant continuum of scattering states. The s.p. Berggren basis is generated by a hnite- 
depth potential, and the many-body states are obtained as the linear combination of Slater 
determinants spanned by resonant and non-resonant s.p. basis states. GSM is a tool par 
excellence for nuclear structure studies which includes all continuum effects and correlations 
between nucleons simultaneously. However, an absence of the projection on the asymptotic 
decay channels does not allow at present for the application of GSM to the description of 
nuclear reactions. 

In this work, we shall formulate the SMEG with the two-particle scattering continuum 
and present the theory of the two-proton (2p) decay. The 2p radioactivity, which has been 


observed in 2002 by Pfiitzner et al. and Giovinazzo et al. is a new spontaneaous 
decay mode in Nature which may appear in even-Z nuclei beyond 2p drip-line. Various 
theoretical models have been proposed to describe this new form of the radioactivity 


The profound relation between the appearance of the 2p radioactivity in even-Z systems 
and the short-ranged pairing correlations gives a hope that in future studies of the proton- 
proton spatial-momentum correlations in the asymptotic state of 2p decay, one will learn 
about the basic features of the pairing held, such as its radial dependence or the multipole 
structure. A possible inhuence of the nucleon-nucleon (NN) correlations in the decaying 
even-Z nucleus on the asymptotic state of two protons, can be extracted from the data 
only in the unihed theoretical framework which contains all ingredients which are necessary 
for a description of the initial state of A-nucleon system, the hnal state of A — 2-nucleon 
system and the asymptotic state of two emitted protons. In the SMEG formalism with the 
two-particle continuum, one can describe: (i) the asymptotic physical two-particle states. 
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(ii) the configuration mixing in the wave functions of the decaying system with A nucleons 
and the daughter systems with A-1 and A —2 nucleons, (hi) the coupling of discrete states in 
the parent system A to the decay channels [(A — 1) 0 (1)] and [(A — 2) 0 (2)], as well as (iv) 
the complete antisymmetrization of the wave-functions in both parent (A) and daughter 
{A — 1) and {A — 2) systems. Hence, the SMEC allows to extract features of the NN 
correlations not only from the width of the decaying state but also from the correlations 
both between the emitted protons and between the emitted proton(s) and the daughter 
nuclei (H — 1), {A — 2). Presently available experimental data allows to address only the 
first of those two theoretical challenges. 

In this work, we shall apply the SMEC formalism for the description of the half-lifes 
for the 2p decays from the ground state (g.s.) of ^®Fe, ^®Ni and ^^Zn. This formalism has 
been applied before for the study of the 2p decay from the excited state I 2 in ^®Ne 0 . 
Here, we shall compare a diproton decay (a direct 2p decay) with a sequential 2p emission 
through either the continuum states correlated by weakly bound states of A-1 nucleus, or 
the resonance(s) of the nucleus A-1. The residual interaction between discrete states and the 
scattering continuum states is the Wigner-Bartlett zero-range interaction what restricts our 
description of the diproton decay with three-body asymptotics to the scenario, introduced 


first in the i?-matrix studies by Barker 2, consisting of the emission of (2p) cluster and 
incorporating the final-state interaction between two protons in terms of the s-wave phase 
shift. 

The paper is organized as follows. The essential elements of the SMEC formalism with 
one particle in the scattering continuum are discussed in chapter H. The extension of the 
SMEC formalism to include both one- and two-particle continuum states will be discussed 
in chapter HI. Chapter IV is devoted to the detailed presentation of the theory of the 2p 
emission. We shall discuss different limiting situations corresponding to (i) the direct 2p 
decay with three-body asymptotics (cf sect. IV.C), (ii) the direct (2p) cluster emission (cf 
sect. IV.B), and (hi) the sequential 2p emission through either the correlated continuum 
in the absence of intermediate resonance (s) or through the resonance (s) in the intermediate 
A-1 nucleus (cf sect. IV.A). All details which are not necessary to follow the main ideas 
of the extended SMEC formalism and the 2p emission theory are put in the appendices. 
Application of the SMEC formalism to the description of the 2p decay from the g.s. of ^^Fe, 
^®Ni and ^'^Zn is discussed in chapter V. The main conclusions are given in chapter VI. 
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II. SMEC WITH ONE PARTICLE IN THE SCATTERING CONTINUUM 

The SMEC describes the nucleus as an open quantum system (OQS) P|. The total 
function space consists of two sets: the space of L^-functions used in the standard SM and 
the function space of the scattering states into which the SM states are embedded. Discrete 
SM states of the A-nucleon system can decay only when their coupling to the scattering 
wave functions is taken into account. 

The two sets of wave functions are defined by solving the Schrodinger equation: 

Hsm^ = ( 1 ) 

for the discrete SM states of the closed quantum system (CQS), and the Schrodinger equa¬ 
tion: 

= 0 ( 2 ) 

d 

for the scattering states of the environment. Here, Hsm is the standard SM Hamiltonian and 
Hcc = Hq + Vcc is the standard Hamiltonian used in the coupled-channel (CC) calculations. 
The channels are determined by the motion of unbound particle relative to the residual 
nucleus with A-1 bound particles in a certain state The states of a daughter 

nucleus are discrete SM states, are scattering states projected on the channel c. The 
channel numbers c are defined by the quantum numbers of the states j of the daughter 
nucleus and those of the unbound particle which are coupled to the total angular momentum 
and parity. The states of the A-1 system are assumed to be stable. 

By means of the two functions sets: Q = V = {|^e)}, one can define the 

projection operators: 

N 

Q = 

i=l 

POO 

P = / dE\^E){^E\ (3) 

Jo 

with 

= 0 

P|4;^> = 0 (4) 
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and the projected Hamiltonians: QHQ = Hqq and PHP = Hpp. We identify Hsm with 
Hqq the CQS Hamiltonian and Hcc with Hpp. The Schrodinger eqnation in the total 
fnnction space is 

{H-E)\'^e) = 0 (5) 

Assnming Q + V = Id, one can determine a third wave fnnction which is the contin- 

nation of the SM wave fnnction into the continnnm V. Fnnction is obtained by 

solving the CC eqnations with the sonrce term ItCj) = HpQ\^i): 

iy+>(J5)>=Gp’(-E)|K>.> (6) 

where 

d;\E) = P{E-Hpp)-^P ( 7 ) 

is the Green’s fnnction for the motion of a single nncleon in the V subspace, E is the total 
energy of the nucleus A and Hpq = PHQ. Using the three function sets: {1$^^)}, {I'Cb)} 
and one constructs the solution |\l/£;) = + -P|fI/£;) in the total function space 

with: 

Q\^e) = {E-nQQ{E))-^HQp\^E) 

P\^e) = \^e)+G^;-\e)HpqQ\^e) ( 8 ) 

One obtains: 

= Ife) + (9) 

i,k 

In the above equations, Hqq{E) stands for the energy dependent effective Hamiltonian: 

HQQiE) = Hqq + HqpG^;\e)Hpq (10) 

in the function space of discrete states which takes into account the modihcation of the 
CQS Hamiltonian Hqq by the coupling to the scattering states. Hqq is therefore the OQS 
Hamiltonian in Q subspace. 

The energies and widths of the resonance states follow from the solutions of the hxed-point 
equations: 

Ei = Ei{E = Ei) 

u = f,(u = Ei) (11) 
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where functions Ei{E) and Ti{E) follow from the eigenvalues of the OQS energy-dependent 
Hamiltonian in Q subspace. The identihcation of E^ and Tj in m with standard spectro¬ 
scopic observables is justihed by an adequate dehnition of the two subspaces Q and P. 

III. SMEC WITH TWO PARTICLES IN THE SCATTERING CONTINUUM 

A. Effective Hamiltonian in Q 

Let us denote by T a subspace of the Hilbert space with the two-particles in the continuum 
and by T the corresponding projection operator. We assume: 

Q + P + T = Id (12) 

which allows to formulate the completness relation in the total function space. Consequently, 
one can decompose the Hamiltonian H into the parts acting in different subspaces and the 
coupling terms between those different subspaces: 

H = {Q + P + f)H{Q + P + f) 

= Hqq + Hqp + Hqj' -|- HpQ + Hpp + Hpp + HpQ -|- Hpp + Hpp (13) 

One can show (cf appendix^ that the effective Hamiltonian T-Cqq{E), which takes into ac¬ 
count the modihcation of the CQS Hamiltonian Hqq by the couplings to the V, T subspaces, 
can be written in the form: 

HQQiE) = Hqq + HqpGUE)Hpq 

+ [Hqt + HqpGp{E)Hpt]G^{E)[Hpq + HppG^{E)Hpq] (14) 

which separates terms due to the coupling with one- and two-particle continuum of scattering 
states. In the above equation, Gp{E) is the Green’s function in T modihed by the coupling 
with V: 

G^{E) = \im[E + ie- Htt - HtpG%{E)Hpt\~^ (15) 

e —>0 

Similarly as in the standard SMEC with one particle in the scattering continuum, we dehne 
the source term: 

\w,) = [Htq + HtpG^;\e)Hpq\\^,) (16) 


7 


and the continuation la;^^) of the SM wave function |<hj) into the continuum T: 

|o;+) = Gi^\E)[HTQ + HtpG^+\e)Hpq]\<!>,) = G^.^\E)\wi) (17) 

Matrix elements of [Hqt + HQpG^i^\E)HpT]G!^\E)[HTQ + HtpG^j^\e)Hpq] in UggiE), 
which correspond to the coupling with T, can be expressed as the overlap between the source 
term (fT?)|l and the function (eq. (ITTIl i. 

One may notice that 'Hqq{E) can be also written in the form: 

nQQ{E) = Hqq + HqtG^^\e)Htq 

+ [Hqp + HQTG^r^\E)HTp\&p\E)[HpQ + HpTG^r^\E)HTQ\ (18) 

which separates explicitely the direct coupling term between Q and T subspaces. Here 
G^p'*{E) stands for the Green’s function in V modihed by the couplings with T; 

G^p\e) = \im[E + ie- Hpp - HptG^p\e)Htp]~^ (19) 

e —>0 

The expression (HHD is obtained by a permutation of projection operators P and T in eq. 

(HI. 

The coupling between Q, V and T subspaces yields also effective Hamiltonians Hpp{E) 
and Htt{E) in V and T subspaces, respectively. For example, the unperturbed Hamiltonian 
Hpp becomes: 

HppiE) = Hpp + HpqGq{E)Hqp 

+ [HpT + HpqGq{E)Hqt]G^p\e)[Htp + HtqGq{E)Hqp] (20) 

in the OQS formalism. Similarly, Htt becomes: 

H.TTiE) = Htt + HtqGq{E)Hqt 

+ [Htp + HtqGq{E)Hqp]G^p\e)[Hpt + HpqGq{E)Hqt] ( 21 ) 

In the following, we shall be interested in applying the SMEC formalism for a particular 
problem of the 2p emission, i.e. the decay from Q to T, for which a relevant operator is 
■HQQiE). 
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IV. DESCRIPTION OF 2P EMISSION IN SMEC 


The effective Hamiltonian Hqq{E) (cf (imi and (ITHll i takes into account couplings between 
Q subspace and the subspaces with one- (P-subspace) and two-nucleons (T-subspace) in 
the scattering continuum. In this expression, all possible emissions of two protons as well 
as one proton are implicitely included. In real 2p decays, certain emission scenarios may be 
less probable, so it is interesting to consider limiting cases of the general emission process 
and isolate appropriate terms describing them. 

If nuclei A and A-1 respect the diproton emission condition considered by Goldansky 

2M- 


Ea-i — ^ 2^^ 


( 22 ) 


then one can suppose that effects of the coupling involving V subspace, i.e. the terms hlgp, 
H^p ... in m and (HHD, are less important. In this case, T-Cqq{E) can be approximated by: 


(S) = ffQQ + ffQTGg>(B)ffTO (23) 

T-Cqq\e) is the effective Hamiltonian in Q subspace describing the direct emission from Q to 
T. If the condition (P|) is not satisfied, then an intermediate system A-1 plays an important 
role in the 2p emission. Neglecting the direct couplings between Q and T subspaces in the 
expression (HI, one obtains: 


= Hqq + HqpG'+'’{E)Hpq 


(24) 


This operator describes emission of two protons through the resonance of an intermediate 
nucleus A-1 and in this case, the emission of the first particle implies automatically the 
emission of a second particle. 

The mechanism of sequential 2p emission process may also occur via the continuum states 
correlated by either weakly bound states of a nucleus A-1, or by weakly unbound states in 
nucleus A-1. Such a physical situation has been studied recently in the 2p decay of the I 2 
state in ^®Ne j^|^. The effective Hamiltonian describing this situation is: 

wg^lB) = Hqq + HqpGH\e)Hpq + [ffgpGG>(B)ffpT]Gi+'(£)[J?TpGG>(£)2/pol 

(25) 
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As before, this expression has been derived from (HHD neglecting the direct couplings between 
Q and T subspaces. The third term in (ESD describes a sequential 2p decay whereas the 
second term is responsible for the Ip decay. In the following, we shall calculate the width 
for 2p emission assuming either a sequential emission, i.e. two succesive and independent 
proton emissions, or a direct emission of (2p) cluster. The intereference between these two 
limiting processes will be neglected in the present studies. 

A. Sequential 2p emission 

The sequential 2p emission may occur either through the resonance of an intermediate 
nucleus A-1 or through the correlated continuum of nucleus A-1. In the following, we shall 
be interested in the latter case. Matrix elements corresponding to the direct Q-T coupling 
are usually much smaller than those of the Q-P coupling due to the smaller Coulomb barrier 
in the latter case. For that reason, we rewrite the effective Hamiltonian (j14jl (cf eq. (j2Sl) as: 

= Hqq + HqpG^^\e)Hpq + HqpGpHptG^^\e)HtpG^^\e)Hpq (26) 

i.e. we neglect Hqp couplings (the direct 2p emission) but keep Hqp couplings (the Ip 
emission). Hence, the interference between the Ip emission and the sequential 2p emission 
decay can be investigated both for close (virtual Q-V excitations) and open (a ’true’ Ip 
decay) Ip emission channels. 

One can see in ()26|1 the source terms for both the Ip emission and the sequential 2p 
emission. Diagonalizing Hqq\e), one obtains energies of states in the nucleus A as well as 
their widths associated with the emission of one and two protons. 

In principle, one cannot separate the partial widths for each of those decay modes. How¬ 
ever, since the couplings corresponding to the 2p emission are in most cases signihcantly 
smaller than those associated with the Ip emission, we shall hrst diagonalize Ti.QQ\E) in the 
SM basis {|<F"^)} neglecting the 2p emission. This provides new basis states {|<l> )} which 
are linear combinations of SM states in Q. Using these new mixed SM states, we calculate 
the 2p emission width for a sequential decay, i.e. we calculate the matrix element: 

6{E) = {^f\HQpGpHpTGi^\E)HTpG^;\E)HpQ\^t) (27) 

describing the sequential 2p emission of the decaying state One should notice, that 
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mixing of SM states due to the sequential 2p emission can be neglected because the 

dominant term is the Ip emission. 

In the following, we shall assume that the subsequent proton emissions are independent, 
i.e. we neglect correlations between the two protons in the continuum and describe the 
interaction of the hrst emitted proton with other A-1 nucleons by a mean-held 
where p denotes the projection operator on the one-particle continuum states. For 
we take a diagonal potential which enters in the CC equations. This implies (cf appendix 

El); 

• Hpp —>• Hq:q: + ph^^^'^^p : Hpp is divided into FTq'q' and ph^^^^'>p. Hqiqi acts in the 
Q' subspace containing the (quasi-)bound states of a nucleus A-1. 

• Hpp —^ Hpipi ph^seq)p . jg divided into Hpipi and ph^^^^'^p. Hpipi acts in the 
subspace V which contains the states of A-1 nucleons out of which A-2 are (quasi- 
)bound and one proton occupies a continuum state. 


• Hpp —>■ Hq'pi ^ Id{A) : i.e. matrix elements of the residual two-body interaction 
involving the hrst emitted proton are neglected. 

With these assumptions, the matrix element S{E) takes a form (cf appendix iBll : 


5{E) = {^f\HQp 


E+ — ph^^'^^p — Hqiqi — HQipi{E~^ — — Hpipi) ^HpiQi 

1 „ 1 




-Hqipi 


X 


E+ - phi^^'^^p - Hpipi 


HpiQi 


E*-phir'p-HQ.Q, 


ffpol-tf) 


(28) 


and the width of a state is given by: 


r(E) = -2Im{6{E)) 


(29) 


The width of a physical resonance state follows then from the solution of the hxed-point 
equations (ED on condition that all subspaces of the Hilbert space involved in the description 
of the resonance decay are dehned adequately. 

In the case of a sequential decay through a resonance in the intermediate nucleus A-1, 
the ehective Hamiltonian is given by the expression (Oil) . With the same assumptions as 
used in the derivation of eq. (I2HD, one obtains (cf appendix ED- 


S(E) = {itlHgp 




— Eqipi{E'I' — phly^'^^p — Hp/p. 


1-1 


H 


P'Q' 


-HpQ^t) 
(30) 
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B. (2p) cluster emission 


In this section, we shall consider the emission of two correlated protons in a form of the 
(2p) cluster. The effective Hamiltonian for this process is given in eq. (I2d|l . Moreover, if one 
includes couplings to the one-nucleon decay channels which are responsible for the external 
mixing of SM states in Q, then 'Hq'q becomes: 

wgj>(S) = ffQQ + HqpG>*\e)Hqp + HQTGfHE)HTQ (31) 


Formally, this expression is derived from eq. m (or eq. m), neglecting couplings between 
V and T subspaces. 


In the following, we assume the two-step scenario for the 2p decay 


221. In the hrst 


step, two protons are emitted in a form of a (2p) cluster. In the second step, the cluster 


desintegrates outside of the nuclear potential of nucleus A-2 due to the hnal state interaction 

(U) of proton 

n 


27|,|28|. The hnal state pp-interaction is taken into account by the density p 
states corresponding to the intrinsic energy U in the proton-proton system 22, |29|. The 


calculation of the density p{U) is based on the s-wave phase shift in pp collision. 
Matrix element of the effective Hamiltonian describing the (2p) cluster emission is: 

S(E) = lM"‘'\HQTG^(E)HT^^t'‘') . 


(32) 


where is the intrinsic state corresponding to the SM state |$i) mixed by either by 

coupling to one-proton decay channels or diproton decay channes (cf sect. IV A 2jl . Working 
with intrinsic states allows to take into account the recoil correction for the daughter 

system A-2. 

Let us consider the completness relation: 




'‘+O0 


dRR^ I drr^ ^ 


\c,r,R){c,r,R\ = R 


(33) 


where c is a channel characterized by the intrinsic state of a daughter nucleus A-2, 
the intrinsic angular momentum of the subsystem formed by two protons R, the spin of 
two protons S, and the relative angular momentum ly between two protons and a daughter 
nucleus A-2. R and ly are coupled to L, L and S are coupled to J 2 p and, hnally, the total 
angular momentum of nucleus A-2 and J 2 p are coupled to J. r is the intrinsic radial variable 
of the cluster and R is the distance between the center of mass of the cluster and the daughter 
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nucleus A-2. Intrinsic state Os of a (2p) cluster is described in the harmonic oscillator basis. 
Since the intrinsic angular momentum of the (2p) cluster is lx = 0, therefore its spin is S' = 0 
due to the antisymmetry of the wave function. 

Sharing of the total energy between the intrinsic energy of cluster and the energy asso¬ 
ciated with the center of mass motion of cluster is taken into account phenomenologically 
by the proton states density p{U), i.e. we suppose that the emission of two protons in Os 
intrinsic state with the intrinsic energy U is distributed according to the density p{U). In 
this case, the problem of 2p emission with the three-body asymptotics reduces to a problem 
with the two-body asymptotics because the degrees of freedom corresponding to the intrinsic 
movement of protons are described phenomenologically by the proton states density p{U). 
The corresponding completness relation (Ea) takes a form: 


(•H-co /•H-oo 

32 / 


dRR^ 


dU 




p(£;)|c(c/), «)«£;), = 


(34) 


c((7)=[t(“‘),0s({/),L,S=0] 

Hence, the matrix element (El can be written as: 


r' + OD /• + CO 

32 / 


dRR'^ 


'0 JO 

/• + 00 /• + O0 


dUp{U){4.t‘>\HQT Y, «){£■(£'). R\GP{E)HTQ\4>t‘^) 

c{U) 


'0 


'0 


dUdRR^wl^iU, R)p{U)J^^J{U, R) (35) 

where Wi^c{U,r) is the projection of the source term \wi) = on the channel c{U): 

R) = R{c{U), R\w) = R{c{U), R\HTQ\^f^') (36) 

The calculation of the source term is given in the appendixO The projected source Wi^dU, 'f) 
does not depend explicitely on U because this dependence in the two-step emission scenario 
follows from the cluster emission process. In the following, the projected cluster source will 
be denoted by Wi^cd)- 

The coupling term Htq is given by the two-body residual interaction: 


Htq = f [YHi] + y] I Q^T\Y q 

\ i i<j / \ i<j / 


(37) 


The Coulomb interaction is included as the average Coulomb held in h and does not enter 
in Htq- 
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The function in (IHK|l is the projection on channel c{U) of the continuation 

= G^rjt\E)HTQ\^f'^*'^) of an intrinsic state in T subspace: 

R) = R{c(U), b (38) 


For the Hamiltonian Htt, we suppose that the total system [A — 2] (g) [(2p)] can be 
considered as a two-body system in the mean-held Uq, i.e. we are considering Htt in a 
following form: 




2/i 


'p(ci) 


(39) 


where jg intrinsic Hamiltonian of the daughter nucleus A-2, and jg the 

intrinsic Hamiltonian of the (2p) cluster. The cluster is described as a particle of mass 
M = 2Mp [Mp denotes a proton mass) and charge Z = 2. P^/2/i is the intrinsic kinetic 
energy of the system [A — 2] ® [(2p)], and /i is the reduced mass of the system. is 

the projection operator on the subspace of cluster states in the continuum of the potential 
F^/2/i -|- Uq. Calculation of the projected function corresponds to solving an 

inhomogeneous Schrodinger equation: 


(£;-4?)|T’> = K> 


(40) 


Projecting eq. (HOD on a channel c{U), one obtains: 

u;l^\u,R) = wUU,R) (41) 

In this representation, pA-2) jg j^pg intrinsic energy of a daughter nucleus and U is the 
intrinsic energy of a (2p) cluster which is distributed according to the density p{U). 


(pA-2) + f/) _ fid) 




C. Direct 2p emission with three-body asymptotics 

In this section, we shall describe direct 2p emission without a simplifying two-step scenario 
in which the hrst step consists of the emission of the (2p) cluster which, subsequently, decays 
due to the hnal-state interaction. Matrix elements of the effective Hamiltonian Hqq\e) (cf 
(j2Sl)) in the SM basis are: 

{4f|H^"j>(i5)|4p = (42) 


14 






where \wi) = Ergl^f) is the source term and = Gip\E)\wj) is the continuation of 

the state in the T subspace. 

Calculation of the energy correction to the SM eigenvalue in (14211 requires a formulation 
for the three-body asymptotic in Jacobi coordinates (see Fig. [T]): 


Xfc \/P'p 


ij 


with /iy 


AiAj 
Ai + Aj 


yk ^(ij)k with f^{ij)k 


(Ai + Aj)Ak 


Ai + Aj -\- Ak 

where Ai = rrii/m, rrii is the mass of a particle i and m is the nucleon mass. The source term 



\wi) and the function are calculated using the expansion in basis of hyperspherical 

harmonics Q . The convenient variables in this basis are: 

• the hyperradius p dehned by : p = ^Jx\ + yl 

• the hyperangle dehned by : arctan (a;fc/|/fc) 

• the ensemble of angles {6xf,,<pxky(^yky^yk) associated with the direction and y^. 

The hyperradius p is the same in all three ensembles of the Jacobi coordinates shown in 
Fig. ^but the hyperangles ak differ in systems Y (the coordinate systems (1) and (2) in Fig. 
^1 and T (the coordinate system (3) in Fig. Q). The hyperspherical harmonics (f^s) 

are dehned as: 


-\/^k’^yk 

xkl 






■Vk i 


a) \^i.k^Xk) 


(43) 
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where ilg = {ak, Xk,yk), and functions (a) are defined using integer order Jacobi 

polynomials: 

(a) = (sin 0 ;)^^'= (cos a)^^''Pn'''‘'''^'^^’^“''''"^'^^(cos(2a;)) (44) 

jg normalization constant, and n = {K — — 4fc)/2- and Yiy^(yk) 

in il are the spherical harmonics associated with x and y, and K is the hypermoment. 
The source term \wi) and the function are calculated in the coordinate system T. To 

simplify notation, vectors (X3, ys) in the dehnition of system T will be denoted by (x,y). 

A decay channel c in the T subspace is dehned as: 


C (t, Kj {J'xi ^y)Li 5, t/2p, t/) 


(45) 


where t is a state of a daughter nucleus A-2, and K is the hypermoment associated with 
the hyperspherical harmonic function y^’i- The angular momenta lx and ly, which are 
associated with the directions x and y, are coupled to L. S is the total spin of two protons. 
L and S are coupled to J 2 p, and the total angular momentum It of a daughter nucleus in 
a state t is coupled with J 2 p to J. The antisymmetrization of the wave functions for two 
protons in the continuum is taken into account by choosing an even value for lx + S. Details 
of the calculation of the source term are given in the appendix O 

The function which is a continuation of the SM state in T, is a solution of the 
inhomogeneous equation: 

{E-Htt)\uj^^'^) = \w,) (46) 


with: 

Htt ~ T 


3>A-1 


p(^-2)+/C + no(A)+no(A-l)+ ^ + V^^\A-1,A) 

i<j 


T (47) 


where is the intrinsic Hamiltonian of the A-2 daughter nucleus, /C is the intrinsic 

kinetic energy of a three-body system: the daughter nucleus and the two protons in the 
continuum. vo{A — 1) and no(A) are the one-body potentials acting on the two protons, 
denoted by an index A-1 and A, respectively. for j > A — 1 corresponds to 

the residual interaction between nucleons of the daughter nucleus and the protons in the 
continuum. — 1,A) -f H*^‘"^(A — 1,A) is the sum of the residual and Coulomb 
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interactions between the two emitted protons. In the basis of hyperspherical harmonics, eq. 
(gni) takes a form of the CC equations (cf appendix lEjl : 




(48) 


where ifcc'(p) is the coupling potential between channels c and d. 

The inhnite range of the Coulomb interaction in does not allow to decouple CC 
equations at inhnity. Consequently, an asymptotic behavior of cannot be dehned 

without an approximation. One way to proceed is to neglect the off-diagonal potentials 
Hcc'ip) for p > On a nd to dehne an effective Sommerfeld parameter from the diagonal 
potentials Hcdp) |3l|- In this approximation, eqs. (HHl) for p > pQ become a system of 
decoupled two-body Coulomb equations. 

If the residual interaction jg ^ contact force: — r 2 ) = I^h(ri — r 2 ) , then 

the contribution to Hcdp) hue to the two-body interaction between two protons in the 
continuum: 


{K, /,, ly, L, J2p, |iC, 4, ly, L, J2p, p) 


cx — Vo / da cos^(q:) sim(a)\k^’^(a) 


Hlyr /(cos(a)/^) 




-3 

p3 


cos a 


)/v^)= 


(49) 


has an ultraviolet divergence for p ~ 0. In this case, solutions of CC equations oscillate with 
a frequency which tends to inhnity as p —0. In general, attractive potentials V{p) ~ p~'^ 
with T > 2 have an inhnite number of bound states 3^ |33|, |3^ and have to be regularized. 
Unfortunately, the standard cutoh procedure for the singular potential at p < po cannot 
be applied since the solution of CC equations would then depend in a random way on the 
value of the cutoh radius po. Hence, the finite-range jg obligatory in solving the most 

general problem of the direct 2p emission. In this case, the CC coupling potentials Hccfip) 
contain a non-local term in p (cf appendix |E)) and, consequently, CC equations (PH|l become 
the integro-diherential equations. Numerical solution of those equations for the problem of 
the 2p decay will be addressed in a future paper. 
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V. DISCUSSION OF THE RESULTS 


Nuclear decays with three fragments in the hnal state are very exotic processes. The 
2p radioactivity is an example of such a process which can occur for even-Z nuclei beyond 
the proton drip line: if the sequential decay is energetically forbidden, a simultaneous 2p 
decay becomes the only possible decay branch. The diproton decay may also be observed 

of the decay 


in a situation where a Ip decay is allowed, as found in the SMEC study 
of I 2 state at 6.15 MeV in ^®Ne j^. However, in this case the diproton decay is strongly 
influenced by an interplay between the external mixing (through the Q — V coupling) and 
the internal mixing (inside of the Q subspace) in SM wave functions [23|, which invalidates 
an idealized picture of an independent decay mode associated with the pairing held. For 
that reason, it is important to search experimentally for the g.s. 2p decay in those nuclei 
beyond the 2p drip-line which are stable with respect to the Ip emission. One should stress 
that the external mixing of SM wave functions is effective also in nuclei with closed Ip decay 
channels, as has been pointed in the studies of the binding energy systematics in sd-shell 
nuclei [l3|. Hence, the many-body states close to the Ip emission threshold can be modihed 
strongly by the residual coupling between Q and V subspaces. Below, we shall study this 
aspect of a diproton decay for ^^Fe and ^®Ni. 

In this chapter, we shall describe spontaneous diproton decays from the g.s. of ^^Fe, 
'^^Ni and ®‘^Zn, which have been observed recently. The valence space used to describe these 
nuclei consists of IsOdO/lp shells for ^^Fe and 0/lp shells for ^®Ni and ^^Zn. As a residual 
interaction between different subspaces we use the Wigner-Bartlett contact force: 


l/l™) = Vo\a + l3F’]S(i^ - ri) 


(50) 


where a -|- /3 = 1, Vb is the strengt 


following, we shall take a = 0.73 and Vq = —900 MeV-fm^ 


1 parameter and is the spin-exchange operator. In the 


A. The decay of ^^Fe 

The 2p radioactivity of "^^Fe has been reported recently [18|, Il9|, |35[ on the basis of accu¬ 
mulated experimental evidence which can be consistently explained assuming an important 
2p decay branch. The reported decay energy is Q 2 p = 1.154(16) MeV. The half-life 
£t of the decay-time spectrum yields ^ T 1/2 = l-S^gJ ms, a somewhat lower value than 
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reported previously 


HQ 


(Ti /2 = ms). T 


le 2 p decay competes with the /d-decay 


and the estimated 2p branching ratio is 0.57(10) The sequential decay through the 
intermediate g.s. of ^"^Mn or the correlated continuum above this state was estimated as less 
based on the model predictions for the Qip which range from -24 keV to -|-10 



1. Calculation of the source term and the function for the diproton decay 


The source term for the direct 2p decay is expanded in the harmonic oscillator basis (cf 
appendix ID| i . The g.s. decay of ^®Fe occurs in the channel: c = Os, L = 0, S' = 0, J'" = 
3/2’''), where is the (intrinsic) = ?)/2^ g.s. of '^^Cr, and L is the relative angular 
momentum between the (2p) cluster and ^^Cr. The internal state of the cluster is Os in the 
harmonic oscillator basis and spins of two protons are coupled to the total spin S = 0. The 
radial wave functions which appear in the source term calculation are gen erated using the 
Woods-Saxon potential with the Blomqvist-Wahlborn parametrization the diffuseness 
parameter a = 0.67 fm, and the radius: Rq = 1.27(74 — 1 )^ 3 . These states are the one-body 
resonances which are regularized using a cut-off procedure j^. The cut-off radius is fixed 
at Rent = 8 fm for lpi /2 , IP 3/2 s.p. resonances and Rent = 7 fm for O/ 5/2 , O/ 7/2 s.p. 
resonances. Those two values of Rent correspond approximately to the top of the potential 
Darrier for those states. The diffuseness of the cut-off function is acut = 1 fm in both cases 

3. 


In the calculation of function we assume that an interaction between the (2p) cluster 

and the daughter nucleus "^^Cr is described by the average potential Uq which is a sum of 
the central Wood-Saxon potential and the Coulomb potential. The parameters of Uq are 
deduced from the deuteron scattering data Q]. The depth of the Woods-Saxon potential in 
Uq is adjusted to obtain the s-wave resonance for a particle of mass 2mp and charge Z = 2 
at the available energy for the 2 p decay. 

The real part of the diproton source function corresponding to the diproton decay of '^^Fe 
(see Table HJ is shown in Fig|21 The imaginary part of the source function, generated by 
an external mixing due to the coupling to Ip decay channels, is very small in comparison 
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FIG. 2; Real part of the diproton source function for the g.s. decay of ^^Fe. 


with the real part. This is due to weak external mixing of different J'^ = 3/2+ SM states 
in ^^Fe. The external mixing is included by considering coupling to the channels ( 
with J'" = 2/ (ground state) and = 3/ (hrst excited state) states of "^^Mn (A-1 system) 
and Ij = Pi /25 P 3 / 2 , / 5/25 / 7/2 waves for the proton in the continuum. 


2. Diproton decay of'^^Fe 


In the calculations for "^®Fe, we take lOKIN interaction (IsOdO/lp shells) in Hqq. 


ion (R 

di y 

ix Conhgurations with the 


This force contains the USD Hamiltonian for the sd shell |42l and the KB interaction for 
the pf shell Q . The cross-shell interaction is the G-matrix 
excitations from IsOd to 0/lp shells are excluded. The results for diproton decay half-lifes 
are summarized in Tabled 


The calculated half-lifes are somewhat longer than found by Dossat et al. 


The 
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Q 2 p (MeV) 

Ti /2 (ms) 

Ti/2 ^ (“®) 

Ti /2 (™s) Qip = —0.1 MeV 

Ti /2 (™s) Qip = 0.05 MeV 

1.138 

21.46 

21.42 

19.80 

19.77 

1.154 

13.33 

13.31 

12.30 

12.28 

1.170 

8.37 

8.35 

7.72 

7.71 


TABLE I: Half-lifes for the diproton decay of the g.s. of ^^Fe for different values of decay energies. 
lOKIN effective SM interaction is used in Q subspace. The strength of the residual interaction 
(ISn|) is Vq = —900 MeV-fm^. The second column corresponds to the approximation of a direct 
2p decay without an external mixing. In the third column, external mixing of SM wave functions 
generated by the Q-T coupling is taken into account in the calculation of the diproton decay. In 
next two columns, we show results including external mixing generated by the Q-V coupling for 
Qip = —O.I MeV (virtual Q-P excitations) and Qip = 0.05 MeV (open Ip decay channel). 

external Q-V mixing of SM wave functions, reduces the diproton half-lifes by ~ 10% for 
chosen values of Qip. Interestingly, this reduction is almost the same if the one-proton 
threshold in '^^Fe is at 100 keV or if the Ip decay channel is opened with an available decay 
energy of 50 keV. One should stress, that this effect depends strongly on the value of Qip (see 
Ref. j23| for the discussion of the external mixing in the decay of I 2 state at 6.15 MeV in 
^^Ne for which Qip = 2.228 MeV) and, hence, its experimental determination is mandatory 
for a full understanding of the 2p decay from the g.s. of ^^Fe. External mixing of SM wave 
functions generated by the Q-T coupling gives a negligible correction to the diproton decay 
half-lifes (cf Table Hj) and can be neglected. 


3. Sequential 2p emission from the ground state of^^Fe 


Half-lifes for the diproton decay of ^^Fe calculated in SMEC using lOKIN interaction 


are longer than reported by Dossat et ah 


35|. One can inquire whether the sequential 


2p emission could play a signihcant role in the g.s. decay of “^^Fe, thereby reducing the 
discrepancy with the data. 

We shall consider the sequential 2p emission through the continuum states correlated by 
the g.s. (J’^ = 2~) and the hrst excited state (J^ = 3“) of “^^Mn. Excitation energy of the 
= 3h state with respect to the g.s. depends on the value of Qip. For Qip = 1.154 MeV 
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in 


44 


Mn, it is E* = 1.455 MeV. In this case, both J'" = 2^ and J'^ = 3^ states are the 
resonances decaying by Ip emission. 

The theoretical scheme for the sequential 2p emission has been described before. In 
a description of the hrst proton emission from the g.s. [J'^ = 3/2^) of ^^Fe we use the 
'45Fe'-reference potential. This potential is generated using the Woods-Saxon potential with 
Blomqvist-Walhborn parametrization the diffuseness a = 0.67 fm and the radius Rq = 
1.27(74 — 1)^/^ fm. The depth Vq of the central part and the strength Vis of the spin-orbit part 
are: Vq = —46.50 MeV and Vg = —8.24 MeV, respectively. This reference Woods-Saxon 
potential is used to calculate the wave functions of s.p. states which are not affected by the 
continuum coupling. For other s.p. states, we have to take into account the correction given 
by the diagonal potential 14c generated by the residual interaction in the CC equations 
To describe emission of the hrst proton to the 2“ and 3“ continuum of ^^Mn, we consider 
Pi/ 2 , P 3/2 and / 5 / 2 , / 7/2 waves for the emitted proton. Let us consider the Ip decay channels: 
(2“,4)^/^'^ with two waves P 3/2 and / 7 / 2 , i.e. cq = (2~, ci = (2“, 47 / 2 )^/^^. The 
diagonal potentials Hcc (c = cq, ci) appearing in the CC equations are renormalized using 


a self-consistent procedure described in j^. We adjust the depth 14c of the central part 
of the Woods-Saxon potential in Hcc for c = Cq, Ci in order to reproduce energies of IP 3/2 
and 0 / 7/2 s.p. states given by the ''^^Fehreference potential. For the remaining channels 
C 2 = ( 2 “,pi/ 2 )^^^^ and C 3 = ( 2 “,/ 5 / 2 )^/^^ involving pi /2 and / 5/2 waves, we take the same 
depths of the Woods -Saxon potential in Hcc as obtained before for P 3/2 and / 7/2 waves, 
respectively, i.e. Wo = V 22 and Vu = Ws- For the channels , we take for 


Q 2 p (MeV) 

Ti /2 (ms) Qip = 0.05 MeV 

Ti /2 (ms) Qip = 0.0 MeV 

Ti /2 (ms) Qip = —0.1 MeV 

1.138 

171.2 

199.6 

258.6 

1.154 

109.9 

127.8 

164.9 

1.170 

71.4 

82.9 

106.6 


TABLE II: Half-lifes for the sequential decay of the g.s. = 3/2/ of ^®Fe for different values of 
the available energy for the 2p decay and different Qip-values. External mixing of 3/2'’' SM wave 
functions generated by the Q-V coupling is included. The strength of the residual interaction 
is Fo = —900 MeV-fm^. 
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Q 2 p (MeV) 

Ti /2 (ms) Qip = 0.05 MeV 

Ti /2 (ms) Qip = 0.0 MeV 

Ti /2 (ms) Qip = —0.1 MeV 

1.138 

235.0 

277.7 

368.9 

1.154 

151.1 

178.1 

235.6 

1.170 

98.4 

115.7 

152.5 


TABLE III: Half-lifes for the sequential decay of the g.s. J'^ = 3/2f of ^®Fe. External mixing 
generated by the Q-V coupling is neglected. Eor other details, see the caption of Table HH 

each wave Ij the same value of the depth of the Woods-Saxon potential in H^c as obtained 
previously in the corresponding channel 

The '^"^Mn'-reference potential used in the description of the emission of the second proton 
from ^^Mn has the same diffuseness as the '^^Fehreference potential and the radius Rq = 
1.27(A — 2)^/^. The strength of the central part and the spin orbit part in ^^Mn'-reference 
potential are: Vq = —47.16 MeV and Vis = —8.37 MeV, respectively. The diagonal potentials 
Hcc are renormalized using the same procedure as described for the emission of the hrst 
proton. The wave functions which are not renormalized by the continuum coupling are 
given by the '"^^Muhreference potential. The one-body operator which appears in 

the expression for the width (cf eq. ()3()|l L is identihed with the Hcc obtained for 
reference potential. 

Half-lifes for the sequential decay, as shown in Tables im and ITm are about one order of 
magnitude longer than half-lifes for the diproton decay and, therefore, cannot explain the 
discrepancy between the reported value j35| and the SMEC results obtained with lOKIN 
interaction. The calculations are performed for Qip = —0.1 MeV (g.s. of ^^Fe is stable with 
respect to the Ip emission), Qip = 0 MeV (g.s. of "^^Fe is at the Ip emission threshold) and 
Qip = 0.05 MeV (Ip emission channel is opened in the g.s. of ^^Fe). In all cases, the 2p 
emission channel is opened. For Qip = —0.1 MeV and Qip = 0 MeV, the sequential 2p decay 
goes exclusively through the ’ghost’ of the g.s. in "^^Mn (the continuum states correlated by 
the proximity of the H = 27 g.s. of ^^Mn) (cf FigE)). The energy centroid of the ’ghost’, 
i. e. the most probable energy of the hrst proton in the sequential 2p decays shown in Tables 
El and EH is < 0 .55 MeV (cf FigEI). This centroid moves slightly with Qip following 

the dependence of the correlated continuum on the position of 27 state in ‘^'^Mn. The full 
width at half maximum of the ghost is = 0.197 MeV, 0.199 MeV and 0.208 MeV 
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FIG. 3: Energy distribution of the first emitted proton in the sequential 2p decay of the g.s. 

= 3/2^ in ^®Fe. The calculations have been performed for different energies of the g.s. J'^ = 2“ 
of ^“^Mn: Qip = —0.1 MeV (solid line), 0.05 MeV (dashed line) and 0.23 MeV (dashed-dotted line). 
The total energy for 2p decay is Q 2 p = 1.154 MeV and the vertical line denotes Qip = Q2p/2. The 
calculations include effect of an external Q-V mixing. 

for Qip = —0.1 MeV, 0.05 MeV and 0.23 MeV, respectively. For Qip = 0.05 MeV, i.e. 
when the g.s. of “^^Mn is inside of the available energy interval [0, Q 2 p] for the seqnential 2p 
decay, the intermediate resonance (J’^ = 2^ g.s. in "^^Mn) contribntion is totally screened 
by the Conlomb barrier and the 2p decay goes throngh ’ghost’ far away from the resonance 
region Qip ± r2-/2. For Qip = 0.23 MeV (the dashed-dotted cnrve in Fig. 121), a fraction of 
the seqnential 2p decay goes throngh the g.s. of ^^Mn decreasing signihcantly half-life with 
respect to the extrapolation of the branch Ti/ 2 (Qip) at Qip < 0.2 MeV. 

The energy distribntion of the hrst proton in the seqnential decay of the g.s. in ^^Fe is 
shown in Fig. |21 for several valnes of Qip. One can clearly see the ’ghost’ of the g.s. in 
^"^Mn at £ ~ 0.55 MeV, which plays an important role in the seqnential 2p decay even if the 
energy of the g.s. in "^^Mn is inside of the interval [0, Q 2 p] (cf the dashed and dashed-dotted 
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curves in Fig. 121 for Qip = 0.05 MeV and 0.23 MeV). In other words, the sequential 2p 
decay for 0 < Qip < 0.2 MeV is predominantly related to the strength of the ’ghost’ at 
e ~ ± and not to the strength of the Ip resonance at Qip. Therefore, the 

sequential 2p decay width in the interval 0 < Qip < 0.2 MeV in the considered example of 


45 


Fe g.s. decay does not reduce to the product of the width for the hrst step (Fip) anc 


branching ratio for the second step, as one would obtain in the semiclassical analysis 


the 


2l|. In 


the interval 0.2 MeV< Qip < 0.275 MeV, the sequential 2p decay the transition through the 
= 2^ g.s. in ^"^Mn is progressively enhanced and becomes dominant for Q > 0.3 MeV. 



FIG. 4: The half-life for the sequential 2p decay from the g.s. = 3/2^ in ^®Fe is shown as 
a function of Qip (circle) along with the half-life for Ip decay (squares). The total energy for 
2p decay is Q 2 p = 1.154 MeV. The calculations include effect of an external Q-V mixing. The 
dashed-dotted line shows the half-life for the diproton decay. 

The ratio between sequential and diproton half-lifes decreases going from negative values 
of Qip to positive ones. External mixing of SM wave functions reduces the half-lifes for 
the sequential decay by about 30% (cf Tables ITTl and Ilf 1|1 . Couplings to the decay channels 
(3“, associated with the excited state ofo^^Mn are relatively important and reduce the 
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sequential 2p emission half-life by ~ 30% for Qip = 0.05 MeV and ~ 38% for Qip = —0.1 
MeV. The above conclusions depend however strongly on the value of Qip in ^^Fe (cf Fig. 

il). 

According to our model, the sequential 2p emission and the diproton emission yield 
comparable half-lifes already for Qip ~ 0.2 MeV. For small positive values of Qip (0 < Qip < 
0.2 MeV) as well as in the case of closed Ip decay channel(s) {Qip < 0), the sequential 2p 
decay half-life changes linearly with Qip (T 1/2 Qip). In this limit, F^p*^^^ S> Fip and the 
2p decay goes via the correlated continuum (the ’ghost’ of the g.s.) in the available decay 
energy window [0,Q2p]- For Qip > 0.3 MeV, one enters in the regime: Fip 3> F^p'^^^ where 
the role of 2^ resonance in the intermediate nucleus ^"^Mn is dominant in the F^p^'^^ In this 
regime, the dependence of T 1/2 on Qip is Gaussian: T 1/2 ~ exp —(Qip — Qi^^^)^- The half-life 
of the sequential emission has a minimum close to Q® ~ Q 2 p /2 and begins to grow again for 
larger Qip-values due to a smaller available energy for the second proton. In the minimum, 
Fip is about 3 orders of magnitude larger than F^p'^'^^ Even though the value Qip ~ 0.2 MeV 
for which a sequential 2p decay becomes comparable with a dip roton decay lies outside of 
the range of Qip-values predicted in various calculations ISll ISSi , nevertheless it is close 
enough to raise doubts about the precise mechanism of the 2p decay from the g.s. of “^^Fe 
which only devoted measurement of masses in this region can remove. 

One should stress that the sequential 2p decay appears always whenever Q 2 p > 0, in¬ 
dependently of the sign of Qip. In other words, the sequential 2p decay competes with the 
diproton mode even for closed Ip emission channels. The separation of 2p decay modes into 
the diproton decay and the sequential 2p emission becomes unjustihed if the two modes yield 
comparable partial decay widths, i.e. F^*^^ ~ T^p^'^^ In this limit, a distinction between 
sequential in time and instanteneous (a direct 2p emission) 2p emission processes looses its 
meaning. Consequently, the hnite-range of residual interactions cannot be neglected and 
the true three-body asymptotics of decaying channels cannot be reduced to any sequence of 
two-body decays. 


B. The decay of ^®Ni 

Recently, a tentative evidence for the 2p radioactivity in ^®Ni has been reported by Dossat 
et ah 0 . In this experiment, one event of 2p decay with the decay energy Q 2 p = 1.35(2) 
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MeV and a partial decay half-life T 1/2 = ms has been identihed j^. Below, we shall 

present the analysis of the sitnation in this nncleus from the SMEC perspective. 


1. Diproton decay of^^Ni 

The g.s. decay of "^^Ni occurs in the channel : c = L = 0, S' = 0, = 0+), 

where is the J'^ = O’*" g.s. of the daughter nucleus '^^Fe, and L is the relative angular 
momentum between the (2p) cluster and the daughter nucleus. Calculation of the source 
term and the function is similar as described above for ^®Fe (see also appendix^. In 
particular, radial wave functions are generated analogously as for ^^Fe and the regularization 
procedure of the one-body resonances is identical as above. 


Q 2 p (MeV) 

ri /2 (ms) (IOKIN)[^, 43, 4^ 

ri /2 (ms) (KB3)[^ 

ri /2 (ms) (GXPFl) 

1.33 

10.3 

11.4 

12.3 

1.35 

6.2 

6.9 

7.4 

1.37 

3.75 

4.2 

4.5 


TABLE IV: Half-lifes for the diproton decay of the g.s. J'^ = 0^ of for different values of 
decay energies and different SM effective interactions. External mixing of SM wave functions is 
neglected. The strength of the residual interaction (1^ is Vb = —900 MeV-fm^. 


Q 2 p (MeV) 

Ti /2 (ms) 

T?/2 ^ (ms) 

Ti /2 (ms) Q\p = —100 keV 

Ti /2 (ms) Qip = 50 keV 

1.33 

10.27 

10.16 

10.09 

10.08 

1.35 

6.18 

6.11 

6.07 

6.06 

1.37 

3.76 

3.72 

3.69 

3.69 


TABLE V: The same as in Tabled but for the diproton decay of the g.s. = 0]^ of ^®Ni. 


For the effective interaction in Hqq we take lOKIN interaction in psdpf shells, as well as 


KB3 


and GXPFl ^ interactions in pf shells. The calculated half-lives (cf Tables HVl 
an dEl) are compatible with the experimental value for all effective SM interactions used 
in the SMEC calculations. The effect of external Q — V mixing on the calculated half-lifes 
(cf Table IV| for lOKIN effective interaction) is negligible for chosen values of Qip. External 
mixing generated by the Q — T coupling is relatively more important than in ^^Fe. 
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2,. Sequential 2p emission from the ground state of^^Ni 

Half-lifes for the sequential 2p decay are shown in Tables IVII and I VI II for lOKIN in¬ 
teraction. We consider the emission through the continuum states correlated by the g.s. 
(J"" = 3/2“) and the first excited state (J^ = 7/2“) of "^'^Co. The residual continuum cou¬ 
pling with Vq = —900 MeV-fm^ reverses the position of 3/2]“ and 7/2/ SM states. This flip, 
which appears for —900 MeV-fm^ < 1^ < —700 MeV-fm^, is a sensitive test of the strength 
of the residual coupling to the continuum states at the beginning of the fp shell. 


Q 2 p (MeV) 

Ti /2 (™s) Qip = 0.05 MeV 

Ti /2 (™s) Qip = 0.0 MeV 

Ti /2 (™s) Qip = —0.1 MeV 

1.33 

25.7 

30.05 

39.7 

1.35 

16.55 

19.3 

25.4 

1.37 

10.8 

12.5 

16.45 


TABLE VI: Half-lifes for the sequential decay of the g.s. W = 0]'' of for different values of the 
available energy for the 2p decay and different Qip-values. External mixing of = 0+ SM wave 
functions via the Q-V coupling is taken into account. The str eng th of the residual interaction 
is Vo = —900 MeV-fm^ and the SM interaction is lOKIN 4^, 43, 4^. 


External mixing reduces the half-life by ~ 15%, i.e. is less important than found in ^^Fe. 
However, the sequential half-life is only a factor ~ 2.5 — 4 times longer than the diproton 
half-life. This ratio decreases going from negative to positive values of Qip. The sequential 
and diproton half-lifes become comparable for Qip ~ 0.15 MeV, in which case the sequential 
and diproton decay modes cannot be considered as independent ones. 


Q 2 p (MeV) 

Ti /2 (™s) Qip = 0.05 MeV 

Ti /2 (™s) Qip = 0.0 MeV 

Ti /2 (™s) Qip = —0.1 MeV 

1.33 

29.4 

34.4 

45.4 

1.35 

18.9 

22.1 

29.0 

1.37 

12.3 

14.3 

18.8 


TABLE VII: Half-lifes for the sequential decay of the g.s. W = O/ of “^^Ni. External mixing is 
neglected. Eor other details, see the caption of Table IVll 















Coupling to the decay channels (7/2“,/^)°^ associated with the excited state of '^^Co is 
between 6% and 9% for all studied cases of the sequential 2p decay of ^®Ni - again signihcantly 
less than found in ^®Fe. 


C. The decay of 


First observation of ^'^Zn and its decay by 2p emission has been reported recently by 


Blank et al. on the basis of both an experimental evidence and consistency arguments. 


The reported decay energy is Q 2 p = 1.48(2) MeV. The experimental partial hal 
2p emission is T 1/2 = ms. The estimated 2p branching ratio is 


.1 

17 


-life for the 


Ml 


The g.s. decay of ^^Zn occurs in the channel : c = Os, L = 0, S' = 0, = O’*"), 

where is the = O’*" g.s. of ^^Ni, and L is the relative angular momentum between the 
(2p) cluster and ®^Ni. The calculation of the source term and the function is similar as 
described above for ^^Fe. 


Q 2 p (MeV) 

Ti /2 (ms) (t=3) 

Ti /2 (ms) (t=4) 

Ti /2 (ms) (t=5) 

1.46 

32.7 

28.4 

27.27 

1.48 

20.6 

17.5 

16.96 

1.50 

12.7 

10.9 

10.65 


TABLE VIII: Half-lifes for the diproton decay of the g.s. = 0]^ of for different values of 2p 
decay energies and different limitations of the Q subspace calculations with the KBS interaction 
The strength of the residual interaction (t^ is Fo = —900 MeV-fm^. 


As an effective interaction in the Q subspace, we take KBS and GXFPl inter¬ 
actions. Results shown in Tables IVIIII fKBS interaction) and IIXI IGXFPI interaction) does 
not include the external mixing, which is very small in this nucleus. Galculated half-lifes for 
the diproton decay are longer than reported by Blank et al. j^. As before, the sequential 
2p decay cannot be excluded a priori. To achieve a better understanding of the 2p decay 
pattern of ®^Zn, one has to know however the Qip values in ®^Zn and ^^Gu. 

Results in Tables Emi and nxi show the convergence of calculated half-lifes with the 
truncation order t which denotes the maximum number of nucleons which are allowed to 
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Q 2 p (MeV) 

Ti /2 (ms) (t=3) 

Ti /2 (ms) (t=4) 

Ti /2 (ms) (t=5) 

1.46 

26.6 

23.0 

22.18 

1.48 

16.3 

14.5 

13.8 

1.50 

10.3 

9.1 

8.67 


TABLE IX: The same as in Table rVTTTl but for the GXPFl interaction 



be excited from the ff /2 orbit to P 3 / 2 , jh /2 and pi /2 orbits, relative to the lowest hlling 


approximation 


We can see that the convergence is nearly attained for t = 5. 


VI. CONCLUSIONS 


All studies of the OQS have shown that the coupling of the system to the environment 
of its decay channels may change the properties of the system P| . These changes cannot be 
neglected when the coupling matrix elements between system and environment are of the 
same order of magnitude as the level distance or larger. In this case, the changes can be 
described neither by perturbation theory nor by introducing statistical assumptions for the 
level distribution. The non-linear effects become important which cause a redistribution of 
the spectroscopic properties of the system. In general, the magnitude of the coupling between 
system and environment depends explicitly on the location of various emission thresholds 
and on the structure of poles of the scattering matrix (S'-matrix). The latter feature is 
absent in the standard SM. For low-/ orbits (/ = 0,1), the coupling to environment becomes 
singular at the particle emission threshold if the the corresponding /-pole of the S'-matrix 
is at the threshold Such couplings may induce the non-perturbative rearrangement of 
many-body wave functions which contain those strongly coupled orbits. For high-/ orbits 
(/ > 2), the coupling is non-singular and can be mocked to certain extent by introducing 
the dependence of the monopole terms in the effective SM interaction on the number of 
particles in valence orbits. Thus, in the same interval of excitation energies, even at low 
level densities, one may hnd coexisting many-body states with largely different susceptibility 
to the coupling to the environment. 

The resonance phenomena are described well by two ingredients. The hrst ingredient is 
the effective Hamiltonian Tigg that contains all the basic structure information involved in 
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the CQS Hamiltonian Hqq. The 'Hqq contains also the coupling matrix elements between 
discrete and continuous states and its (complex) eigenvalues determine both the positions 
of the resonance states and their widths. The second ingredient is the unitarity of the 
S'-matrix which causes a non-trivial energy dependence of the coupling matrix elements 
between resonance states and continuum. 

The SMEC, which has all those features, has been extended in this work to describe the 
decay of the CQS due to the coupling to the environment of decay channels with one- and 
two-nucleons in the scattering continuum. This new theory has been applied here for the 
study of the 2p decay from the g.s. of "^^Fe, ^®Ni and ^'^Zn. The results are very sensitive 
to the Q-values which for these nuclei are either not known experimentally with a sufficient 
precision (Q 2 p-values), or are unknown as in the case of Qip-values. Keeping this warning 
in mind, one finds that calculated values of the diproton width in these nuclei agree rather 
well with the experimental predictions, in particular for even-even ^®Ni and ^^Zn nuclei. 
Different effective SM interactions (cf Tables ITVl IVTTTI and ITXjl give a similar results to 
within ~ 10 — 20%. External mixing, generated by the Q-V coupling, modifies the diproton 
decay width by < 10% in “^^Fe and ~ 2% in ^®Ni^e. the diproton mode seems to be rather 


close to a ’pure’ mode predicted by Goldansky 


24| . This is in contrast to conclusions from 


the analysis of the diproton decay of I 2 state in ^®Ne in Ref. j23|. One should mention, 
however, that the separate treatment of sequential and diproton modes becomes invalid 
when — ^ 2 *^^- In this case, the 2p decay has to be described using a full Hamiltonian 
m with the three-body asymptotic. In the studied case of ^®Fe and ^®Ni, this happens for 
Qip ~ 0.2 MeV. The Q-T coupling provides a negligible contribution to the external mixing 
even in the region of Qip values where Tip -C r^p*^^ 

Somewhat worse agreement between calculated and reported diproton widths is found in 

his difference 
. One should 

stress again that this conclusion depends strongly on the assumed values of both Q 2 p and 
Qip. For Qip > 0.22 MeV, the half-life for a sequential process becomes shorter than the 
half-life of a pure diproton emission. 

Recently, R-matrix approach has been applied for a description of the diproton decay 


“^^Fe. Neither external mixing nor the sequential 2p emission process explain 
if the values of Qip are taken close to the limits suggested in Refs. 1^ 1^ 


in ^^Fe 


221, “^^Ni and ®^Zn In this model, external mixing is neglected and the 


emission process is described by a simple R-matrix ansatz. Smaller importance of the 


31 







external mixing in those nnclei leads to a good agreement with SMEC resnlts. For GXFPl 
effective interaction, one Ends: ^ 7.4^21 and ^ ms in ^®Ni, 

and = 13.8^1^ ms and = io:^7 542n. 

In our studies, we have employed a contact force with the the spin-exchange for the 
residual interaction between Q, V and T subspaces. As a consequence, the three-body hnal 
state in the direct 2p decay could not be calculated (cf sect. II V Gfl and we used the two- 
step scenario to describe this decay. Application of the hnite-range residual interaction 
allows in future to describe this most general case using the theoretical formalism which has 
been presented in this work. 
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APPENDIX A: EFFECTIVE HAMILTONIAN IN Q SUBSPACE 

Solution I'k) of the Schrodinger equation in Q -|- T* -f T = Id, has components in each of 
the considered subspaces, i.e. |\k) = I^/q) -|- |\kp) -|- [H't)- Hence, the Schrodinger equation 
can be written as: 



(E-H) [|Tq) + |Tp)-f |Tt)] = 0 

(Al) 

Projecting (lAlIl 

on Q, V and T, one obtains: 



{E — Hqq)\'^q) = Hqp]'^ p) + HqtI'^t) 

(A2) 


{E — Hpp)\^p) = HpqI^q) + i/pT Tp) 

(A3) 


(E — Htt)\^t) = HtqI^q) + i/pp Tp) 

(A4) 

The component 

4/p) can be found from (jA3jl : 



4/p) = Gp{E) [ifpQ 4 /q) -|- i7pp 4/p)] 

(AS) 


where Gp{E) is the Green’s function in V (cf eq. ([7j)). Hence, Ithp) can be expressed in the 
form (cf ()A4|1 . (IA5|l b 

I^t) = G+(F) [HtqI^q) + HtpGUE)Hpq\^q)] (A6) 
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where Gi^{E) is the Green’s function in T, modihed by the coupling with V (cf eq. 

Using eqs. (HI, (Hi and (Hi, one can demonstrate that I^^q) is a solution of the equation: 


{E - Hqq)\^q) — HqtG^{E) [Htq + HtpG^{E)Hpq\ 


+ HqpG+p{E) 


Hpq + HptGUE) [Htq + HtpGUE)Hpq] 


I^q) 


The above equation can be rewritten in a form: 


(E-Hqq{E))\1>q}=0 


(A7) 


where 'Hqq{E) is the energy-dependent effective Hamiltonian in Q\ 

Hqq(E) = Hqq + HqpG*{E)Hpq 

+ [7/qt + Eqp GUE)Hpt] G+{E) [ Htq + HTpGp{E)HpQj 


(A8) 


APPENDIX B: MATRIX ELEMENTS FOR THE SEQUENTIAL 2P EMISSION 


A part of the effective Hamiltonian Hqq^ in Q which describes the sequential 2p emission 
is (cf sect. pV Afl h 

1 


H 


QP- 


-HpTGp{E)HTpGp{E)HpQ 


(Bl) 


E+ - Hpp - HptG:^{E)Htp 
We assume that the hrst emitted proton interacts with the remaining A-1 nucleons through 
a mean-held This implies: 


Hpt = P 


Htp = T 




A 


J<j=l 


f^Q' 


P' 


' A-l 

■ A-l 


P' ® Id{A) = Hq,p, ® Id{A) 

(B2) 

Q' ® Id{A) = HpiQi (g) Id{A) 

(B3) 


where Q' and P' are projection operators on the subspaces Q' and V' of (A-l)-nucleon states. 
Id{A) is the identity operator which acts on the hrst emitted proton. In Q', all A-l nucleons 
are in (quasi-)bound orbits. In V', A-2 nucleons are in (quasi-)bound states and one proton 
occupies a continuum state. Hence, after the hrst proton emission, all states of the system 
are separated in the two subspaces, denoted P'®p and Q'(g)p, where p is the space of s.p. 
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states in the continuum. The particle occupying the state in p interacts with particles in Q 
and V by the one-body potential Hence, we have: 

A A 

^h(z)+ 


Hpp = P 


= Q' 


i=l 

'A-1 


i<j=l 

A-1 


J2h{i)+ E ''"'"’('.i) 

i<j=l 


i=l 


Q' + ^ ^ p^iseq)p (34) 


where p is the projection operator on the space p. Similarly, one obtains: 

HtT = Hp:p: + 


(B5) 


Hence, the expression (IBljl can be written as: 




1 


QP- 


E+ — ph^^^l^p — Hqiqi — Hqipi [£'+ — — Hpipi] Hpiqi 

^ ^Q'P'-PZ - PTEPPT, -E- ^P'Q'-PZ -ZTTEEiT-E- ^pQ 


= Hr 


E+ - - H 

-Ho'i 


p,p, E+ - ph(^^i)p - 

-H, 


Q'Q' 


E+ - E+ - pE^^<i)p - E+ -ph(^^<i)p - Hq>q>^^^ 

M in the above equation stands for: 


M = Hq,q,+ph^^^^'^p + Hq,p,- 


-Hpiqi 


(B6) 


(B7) 


EP — ph^^'^i^p — Hpi pi 

Let us now calculate the matrix element of the operator in ()B6fl for the state Let 

us dehne the completness relation in V: 


P+OO 


E/ de\t,e,l,j-,J}{t,e,l,j-,J\^ li 


(B8) 




where \t) is an eigenvector of Hqiqi corresponding to an eigenvalue Et . The state \e,l,j) is 
an eigenvector of ph^^'^^'^p with energy e (e > 0 ), the orbital angular momentum I and the 
total angular momentum j. Discretization of the energy integral in (IB8fl yields: 

ri,=+oo 

^ ^ ^ ^ I^7 (^7 I 

c n=l 

where the channel is dehned by c = {t, l,j] J). In the above expression, Ag is the discretiza¬ 
tion step and the states |e„,/, j) are normalized as follows: 

^n. n.l 


{^ni^i j\(^n' I^i j) 


-'n,n' 

a: 
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Inserting four times the completness relation ()B9fl in (IB6fl . one finds: 


+ 00 




ci,C2,C3,C4 ni,n2,n3,n4=l 

X {C2, Cnal-ffQ'P' 


E~^ — — Hpi pi 

x{c4,e„JhfpQ|<lf)Ag 


^^'E+ - M 

HpiQi\c 3 

5 )(C3, ^n3 ' 


^2) Cn2 ) 
1 


E+ - pE^^l^p - Hqiqi 


1^45 6714 ) 

(BIO) 


Let us first consider the matrix element: (ci, \E~^ — M\c 2 , Cnj), which can be written as: 


{ci,eni\E+ - M|c2,e„2) 

= {h, J\E+ - ph^^^^^p - Hqiqi - HQipi 

{E Ef^ 6ni )<^ci,C2*^eni 


E+ - ph(-^^i^p - Hpi 


-HpiQi\t2, e„2, l2,j2] J) 


,6712 A 

1 


E+ - - H. 


-HpiQI 1 ^ 2 ) 5(6771 ,h,jl),ien2 Hj2) 


ni — ripipi 

Here, |fi) is an eigenvector of Hqiqi with an eigenvalue Eq and: 

^ r +00 


A, 


(Bll) 


-I _ r-roo 

{ti\HQipi— - - - HpiQi\t2) = / drwty,ppj'{r)ujl'^J^^,~y"^\r) (B12) 

E+ - e„i - Hpipi Jq 


t', V and j' in jnn denote the bound state of a daughter nucleus, the orbital angular momen¬ 
tum and the total angular momentum of the second emitted proton, respectively. 
denotes the projection of the source term on the emission channel: d = {f, JtJ, where 
Jq is the total angular momentum of |fi). is the projection on the same channel 

of the continuation of the state 1 ^ 2 ) in the continuum. The expression ()B12f) corresponds to 
the emission of a second proton with energy E — from a nucleus A-1. For > E, the 
second emission is impossible and, consequently, the expression (EH becomes real. The 
source term and the continuation of the SM wave function in the continuum is calculated 
similarly as in the standard SMEC. Thus, the matrix element (IBllfl takes a form: 


{ci,e„i|E+ - M|c2,en2) 
(L' ~ L'ti ~ 6ni)5ci,c2 5e,7i 


,cn2 


E 


f + oo 


""1^)5(6771,/!^ 


1 


The inverse operator M) in the matrix element (cf eq. (Eini)) is obtained 

by diagonalizing E~^ — M for each energy 6^ in the basis of Q'. 
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The term ( 03 , 6 ^ 3 |(E+ Hq^q/) ^|c 4 ,e„J in (IBlOp equals: 

{c3,en3\{E+ - - HQ>Q,)-^\Ci,en^) = [E - - ^i3]“^^(c3,C4)^(e„3,e„J^ (B13) 


and the term (<h^|ifQp|ci, e„j) is: 


{^f\HQp\ci,en,) = / dr 


(B14) 


Here, i ^ (^) is the source term projection for the emission of the first proton 

and is the radial wave function of the state |e„j,/i,ji). The calculation of 

(c 2 ,e„ 2 |i/Q/p/(E+ - - Hp,pi)-^Hp /Q/ 103 , 6 ^ 3 ) is identical to the calculation of the 

term in eq. (inni). Hence, the matrix element becomes: 

H-CXD /• + 00 

rJr 


= E E / 

*11*2 ,*3 iti 11=1 


dr wlA.^ ij{r)ue„,i,j{r) 


X {t 




Hq'Q' — Hq:pi[E+ — Cn — Hpipi] ^Hpi 


Q’ 


-\t2) 


r +00 

X / 


'•+CXD 




E+ -e. 


n 


(BI 6 ) 


o u 

and the partial width for the sequential 2 p emission is given by eq. (I 2 ni). 

For e„ > E, the emission of a second proton is impossible. Hence, the contribution of 
different terms with is real. Since we are interested in calculating the emission width, 

therefore the energy summation in ()B15fl is restricted to an interval from 0 to E. 

If the sequential decay occurs through a resonance in the A-1 intermediate nucleus, then 
one has to consider the operator HqpG^p\e)Hpq (cf eq. fl^ b Assuming (jB2j) - (jB5j) . one 
obtains: 

1 


H, 


QP- 


-H 


PQ 


E+ - ph^seq)p _ [E+ - ph(^^^)p - Hp,p,] ^ Hp.Q, 

Inserting two completness relations (cf eq. ()B9fl h one finds the matrix element: 

n-\-oo 


(BI 6 ) 


n +00 

*1,*2,*,7 n=l 


X {t 


X 


E^ — Cn — Hqiqi — iFQ/p/[F^+ — Cn — Hpipi] ^Hpiqi 


1 ^ 2 ) 


r‘+oo 


(B17) 


from which the partial decay width can be calculated. 
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APPENDIX C: THE SOURCE TERM FOR THE EMISSION OF (2P) CLUSTER 


Projected source term for a direct emission of two protons as a cluster is: 

Os, {Lrel, S), J2p; J, R\HTQ\cl>t'^) (Cl) 

where c = Os, {Lrei, S)', J 2 p] J) is the decay channel in T. is the intrinsic state of a 

daughter nucleus. Os is the intrinsic state of a cluster, L^ei is the relative angular momentum 
between the cluster and the nucleus A- 2 , and S is the spin of the cluster. Lrei and S 
are coupled to J 2 p, and J = J^{m) + J 2 p is the total angular momentum of total system 
[A — 2] ® [2], R in (IClf) denotes the relative coordinate between the daughter nucleus and 
the cluster. The source term Wi^dR) is localized and can be developed in the harmonic 
oscillator basis: 


Wi^dR) = R ^rel^rel ]J2p]J] Nrel\HTQ\(d'^d 


Nr 


where TZjsf^^^L^^dR) i® harmonic oscillator wave function characterized by the radial quan¬ 
tum number Nrei and the angular momentum Lrei- 

In the formalism of second quantization, the coupling operator between Q and T (cf eq. 
(|nZl)) is: 


htq - - 

a < p 
'Y < (5 
r 


a/1 + y/1 + '^7,(5 




(C2) 


where ^ 5 is the antisymmetrized, reduced matrix element of the residual interaction. 
a, ...6 are the eigenstates of the one-body potential with an origin in the laboratory frame. 
Since Htq is expressed in the laboratory frame and Wi^dR) is calculated in the frame asso¬ 
ciated with the relative coordinate i?, therefore one has to change the coordinate system in 
order to calculate the projection of the source. Let us consider the non-spurious SM states 
\ti) and which correspond to intrinsic states and 10^*/. By dehnition, these 

states can be written as: 


l*i) = lAblKot") . = (C3) 

where and |‘hoo) ground states of the center of mass of nuclei A -2 and A, 

respectively. Let us consider the following matrix element expressed in the laboratory frame: 

{HQsdL,S)-hp.dN\HTQ\^t) (C4) 
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N and L in (IDl are the oscillator quantum numbers characterizing the state of a cluster 
with respect to the origin hxed in the laboratory frame. Using Moshinsky transformation, 
one can write (El as follows: 

(OOiVL; L\Ng, Lg, N^el, Uef, L){Ng, Lg, Uel, s, J, Nrel\HTQ\<!>t) 

Na,La,N'rel,Lrel 

where {OONL; L\Ng, Lg, Nrei, Lrei; L) is a Moshinsky coefficient. Ng and Lg in the above 
equation are the quantum numbers corresponding to the motion of total system with re¬ 
spect to the laboratory frame. The residual interaction does not act on the center of mass 
coordinates of the system, what implies: 

Ng = 0, Lg = 0, N^ei = N, L^ei = L . 

One obtains: 


(L, Os, (L, L); J 2 p; J, N\HTQ\^f) = (OOiVL; L|00iVL; (L, L); hp, J, N\Htq\4"'^) 

(C5) 

Hence, the matrix element in the laboratory frame has been transformed into the matrix 
element in the relative coordinates. Using the analytical expressions for the Moshinsky 
coefficients, one can rewrite CD as: 


w. 


. N {2N+L)/2 

= E ( 1 Os, (L, L); hp, J, N\HTQ\^f) 


N 


(C6) 


where: Unl{R) = RRnl{,R)- Using the Wigner-Eckart theorem, one writes Wi^dR) as: 

1 f A ^ {‘2N+L)/2 


w, 


AR) = - E 


(q < /3) G cont 
(•7 < 5) ^ disc 
r, N 


\/l -|- ^ 7 ,<5 ^ ^ ^ 


UNliR) 


xJ 


Jti J 2 P J 

J 0 J 

r r 0 


11 [a,Ss|'’| |4p (Os, L, S,J 2 „N\\ (afif) ((a, /3)^| V| | ( 75 )') (C7) 


Let us consider a term: 


E 


{(y<P)£cont 




: (Os, L, L, J 2 p, N\ I {ai3f){{a,f5f\\V\\{j6f) 
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in eq. dCZl). Inserting twice the completness relation defined with the harmonic oscillator 
wave functions and integrating over the energy of states a and (3, one obtains: 


^ f 2(0s, L, S, Jap, iV| (ni, hji, ria, k, ^ 2 )'") ((^ 3 , h, is, h, ji)^\V\{'ySf) 

(n-l, < (tI 2 , i 2 . J 2 ) 

(ns, 13,33) < (" 4 , U, 34 ) 

1 1 

X - qi^jjj3){j2\l - - qi^,jjji){j2\l - qipjfslh)} (C8) 


where ji represents the state and similarly for j2,j3,j4,- qia,ja qi^j^ are the 

projectors on proton (quasi-)bound states having quantum numbers la-ija and respec¬ 

tively. The term (Os, L, S', Jap, iV| (ni,/i, ji, na, ^2, ^2)^) in equals: 


^2(1 + 


.kkSLil - (-1)^+1+') < 


31,J2) 


h Si ji 
h Sa 32 
L S T 


> (ni, /i, na, /a, L\N, L, 0, 0, L) (C9) 


The term {{n 3 ,l 3 , j 3 ,n 4 ,l 4 , ji)^\V\{'jS)^) in ()(18|1 equals: 


i<j 




1 

iy/2(l -h (5ij) y/ 2(1 5j3j4) 


(1 - - (-1)'52+^2+1) 


L 2 S 2 "^^2 ^ 2 ^ 2^2 





< > 


ijS Si Jjg 


'^1 ji 

X < 

^4 jji 


h '^2 jj ^ 


Li Li T 

V. J 


L 2 S 2 r 


*■ jijjL2S2 


^ ’ ^i4! ^i4! -^11-^15 -^1) , Ti) X Tij , Ij , Ta |.^a^ *^2’ ^ 2 > ^ 2 > -^ 2 ) 

X (iV(, L;, <, Li, Si, r|f"|iV', L', n'a, Ta, Sa, T) 


(CIO) 


Since the residual two-body interaction does not act on the center of mass of two 

particles, we have: 

7V( = iV' , L[ = L' . 

Moreover, in the case of Wigner-Bartlett interaction which is used in this work, the spin 
(S'! = Sa) and intrinsic angular momentum (/( = I 2 ) are conserved, what implies: Li = La- 
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Hence, takes a form: 

I ( 7 , Sf){{j3, J4f\v\ {i, jf) 

i<j 


i<j Li,Si,N[,L[,n[,l[,n'^ 




J3 J4 


^is jjs 
hi ^‘2 jji 

Li Si r 




X 


^ jijj^iSiinj^lj^^ hit -^i|-^i) ^it ^it ht -^1) 


h "Sl ji 
h ^2 jj 
Li Si r 

X {n,k, n,, Li |iV{, L',, n', Li) (iV(, L;, n'l, l[, h, l[, Li, Si) 


(Cll) 


Inserting (j(^8jl . and (Unu in eq. (IHtI) . one obtains: 


Wi,c{R) = - ^ unl{R) 

T ^{')<S)^disc^N 


A 


(2N+L)/2 


X 


E 


1 


A-2^ 

1 1 


j(rf 


Jti J2p J 

J(fyA 0 J 

^ i 

r r 0 


:(ti||[a>~5]^||$f) 


1 + ,79 1 + 5 


J1J2 ^ 


7,5 


(la,3a) < (Ip,op) 
(911, ll, 3 l) < (" 2 , 12,3-2) 
(ns, 13,33) < (" 4 ,4,44) 


X {(ji|l - qi^,jMz){j2\l - qipj^ji) - (-1)^ - g«„,i„b4)(j2|l - qi^jp\h)] 

h Si ji 
h S 2 32 


xjij2LSV2 


> (rii, /i, n2, k, L\NL00, L) 


L S T 


xE 

*<j 


=2(ij|(7,5)^) < 


^i3 jjs 

hi "^2 jji ^ +3+4-^1*^1 ^ 

Li Si r 

/ 

xj+j+i5^i{n 

i3 ha t +'4 1 hi t LijNi,Lhn[,h,Li} 




'5l + 

h "^2 + 
+1 5i r 
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APPENDIX D: THE SOURCE TERM FOR THE DIRECT 2P EMISSION 


The projected source term for the direct 2p emission is: 

WiM = K, (/,, ly), L, S; J2p; J, p\HTQ\<!>f) 


(Dl) 


where the emission channel is defined as: c = (t, K, ly), L, S; J 2 p; J) (cf sect. II V Cfl . Using 
the second quantization form (ina for Htq, one obtains: 


WiM = K, '»). S; J2p; J, p\ Y. Uai,; 




a < 0 
■7 < 5 

r 


V^l + \/l + S-t,5 


Since the state represented in bra has two protons in the continuum and in the ket all 
nucleons are in (quasi-)bound states, the operators can only annihilate particles in the 
continuum states which are present in the bra vector. Hence: 


Wi,ciP) 


-p^/^{t,K,{l,,ly),L,S-,J2p-,J,p\ 

ct < 0 G cont 
7 < 5 G disc 

r 


yi+ 5ij,,j3 yi+<5^,5 




Separating the part of the operator acting on continuum states from the part acting on 
(quasi-)bound states and applying Wigner-Eckart theorem, one Ends: 


WiAp) = -p 


— _n5/2 




--J 


a < 0 G cont 
7 < (5 ^ disc 

r 


It J 2 P J 
J 0 J 
r r 0 


> (t||[a>~5]^||<hf) 


a/I + <^7,5 

X {K, ly), L, S; hp, p \I (a, 0)^) ((a, 0)^\\V\\ ( 7 , 5 )^) 


(D2) 


Let us consider the sum over energies of proton states a and 0 in the above expression: 

Y {A', L, S; J2p, p| |(a, ; 3 )'’) {(a, ; 3 )'’| I v^l I (7, Sf) 

a<f3Gcont 

= EE {K, {h, ly), L, S; J2p, p\\{a, 0)^) {{a, 0)^\\V\\{^, 6)^) 

a<p ea,ef3 

= (4, ly),L, S; J2p, p\ I(U - fe) ® (1/3 - qp )I |l"l 1(7, (D3) 

a<P 

where a denotes the quantum numbers la, ja and 0 the quantum numbers Ip, jp. The 
operator is the projector on proton (quasi-)bound s.p. states with the orbital angular 
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momentum and total angular momentum ja- In a similar way, projects on the subspace 
of proton(quasi-)bound states with the angular momentum and total angular momentum 
jj 3 . Hence, the projected source term Wi^dp) can be written as: 


Wi,c{p) = -P^^^ 


--J 


OL < 0 

7 < <5 G disc 

r 


It I^2p J 
J 0 J 

r r 0 




\/l + ^7,(5 

h h u 

x{K, {l^,ly),L,S; J2p, p\\{la - Qa) ® ( 1 ^ - 11^11 ( 7 , 

The term (1 q — g^) ® (1^ — q^) in (jD4jl can be written in a following form: 

(la - Qa) ® (1^ - qp) = la ® Ig - la ® gg - ® ® 


(D4) 


(D5) 


The sum over quantum numbers a and (3 yields for the hrst term in (insi) the identity in the 
space of two-particle states: 


la (g) Ig = 1 (g) 1 

a</3 

Hence, the contribution of this term in Wi^dp) (eq. (ID4ll j can be written as: 


(D6) 


7<(5srfisc;r 


J 


Pt J2p J 

J 0 J 

r r 0 


W(| I [a,aj]‘’| |4>P(A', (i„ l,),L, S; pWVWi-,, Sf) 


(D7) 


States 7 and 5 are (quasi-)bound and can be expanded in the harmonic oscillator basis: 

= E {nd^j^nshjs; r| ( 7 ^)^) \n-y, d, j^ns, I5, js] T) (D8) 

nj<n§ 

where and ns are the radial quantum numbers associated with the harmonic oscillator 
expansion. Passing from j-j to L-S coupling scheme and applying Moshinsky transformation, 
one rewrites expression ()D8|I in a following form: 

1(7^)^) = E (’^7^7^445; ri( 7 ( 5 )^) y^ 7 A',s'(’^ 7 >^< 5 ) 

n^<ns N',L',n',l' X',S' 

X {nd^nsls\N'L'n'l'-,X')\N'L'n'l'-,X'S') (D9) 
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where the coefficient is given by: 


Vy,S'{n^,ns) = 


1 + (-1)'^'+^' 


X j-yjsS'X' < 


Iry S-y 

h ss js 
X' S' r 


Introducing the expansion ()D9jl in the following term of (|D7|) ' one obtains: 


(DIO) 


{K, (/,,g,L,^;J2p,p||D||(7<5)^) 

nj<ns N',L',n',l',\',S' 

{n-yl^nsl 5 \N'L'nl'-,X) x {K, L, S; J 2 p, p\\V\\N'L'n'l'] X'] S'-,T) (Dll) 


where matrix elements of the residual interaction are: 


{K, (/,, /,), L, 5; J2p, p| |D| liV'LV/'; A'; 5'; h) 

= f / dadMf^5Cos2(a)sin2(«)vh^'=’'^'=(a)y^:(4)xM 

v^psin(a) 


x'JZ 


N'U 


y/Th 


T^n'V 


pcosya] 

'\/‘^Px 


Yv{^x)x'si^s)^ly,L' 


pcos(a) 

y/l^x 



(D12) 


Px and Py in ()D12fl are the dimensionless reduced masses associated with different sub¬ 
systems related to x and y directions in the Jacobi coordinates system T (cf sect. IIV (Jl . 
J^v'L'(v^psin(a)/y//i^) and TZn'v (p cos(q;)/ y/2px) are the radial functions in the harmonic 
oscillator basis associated with the center of mass motion of two protons with respect to the 
nucleus A-2 and the relative motion of those two protons, respectively. 

Contribution of the second term in ()D5fl to tCi,c(p) is: 


E 

a < /3 

7 < 5 ^ disc 

r 


's/I -|- '^7,1 


--J 


It J2p J 

J 0 J 

r r 0 




X (g (/., /,), L, J 2 „ p\\U®q-p\\V\\{l, 5f) 


(D13) 


The last term in (jP 1 3jl is calculated by inserting two completness relations in the harmonic 
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oscillator basis. One gets: 


{K, {l,,ly),L,S-J2r>^p\\l^®rp\\V\\{l^5f) 


= E E (■^5 ^y)l -^5 P\\ni,h,jin2,l2,j2) 

h^Jl n[, 

^2 » ^2 ? J2 ^2 ’ ^2 ’ ^2 


x(ni,li,jin2,l2,j2l 


la 


all^l,^i,J>2,^202) X K,^l0>2,^202ll^ll(7,<^)’") 


(D14) 


To calculate the overlap (K, (I^, ly), L, S; J 2 p, pjjni, Ii,jin 2 ,12,j2}, we make first a transfor¬ 
mation from j-j to L-S coupling scheme and then use the Moshinsky transformation for a 
state Ini,li,jin 2 ,l 2 ,j 2 } (cf (1091)). 

Contribution of the third term in (jnSI) to tCi,c(p) is: 


E 

6l < (3 

7 < 5 G disc 

r 


^1 -\- 


--J 


It J2p J 

J 0 J 

r r 0 


> (t||[a>~5]^||<hf) 


X {K, (4, ly), L, J 2 p, p\\qa®lp\\V\\ ( 7 , 5)^) 


(D15) 


This term is calculated in the same way as the term (Iml . 

Finally, the contribution of the fourth term in (ID5jl to tCi,c(p) is: 


E 

6l < (3 

7 < 5 E disc 

r 


\/l + ^7,< 


--J { 


It J 2 P J 

J 0 J 

r r 0 




X {K, (/., ly), L, S; J2p, pm®qp\\V\\ ( 7 , 5)^) 


(D16) 


To calculate this term, we insert two completness relations in the harmonic oscillator basis. 


APPENDIX E: CALCULATION OF THE FUNCTION ujM FOR THE DIRECT 
2P EMISSION 

The state which is the continuation of in T subspace, is the solution of the 

inhomogeneous Schrodinger equation (gni)- can be expanded in the channel represen- 
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(El) 


tation (in the Jacobi coordinates system T) as: 

C 

= [ dp\c,p)ujl'^J{p) 

C 

where 00 j^\p) is dehned as: 

ui‘j*J{p) = p‘'Hc,p\u,‘j*>} (E2) 

Inserting JEB into the eqnation (gni), and then projecting on the channel c, one hnds: 

pHc, p\E - Htt\Y^ f dp'p'^c'p')uj\'^J{p') =Wj,cip) (E3) 

c' d 

where Wj^dp) is the projection of the source term Wj on the channel c (cf appendix |D)) and 
Htt is given in (gZl). By dehnition, belongs to T. Hence, one can rewrite ()E3f) as: 

pHc,p\f{E ‘ip'p'’|c'p')dtV') = wM (E4) 

d 

The projection operator T assures that the particles in the continuum do not occupy the 
(quasi-)bound states of the daughter nucleus. 

In the following discussion of this appendix, we shall leave out T and replace it by the 
identity operator I^. The effect of the projection operator T will be taken into account 
subsequently by the method described in the appendix |Fl Let us dehne by the 

solution of equation ()E4fl without the projection operator T . One obtains: 


(c, p\{E-H)\^ f dp'p'^^\cp)d^]’'^{p) = Wjdp) 


(E5) 


where d^^P{p) is: 


‘^Pc’“(p)= P^'‘‘{‘:<Ppp^'°) 


(E6) 


1. Calculation of the term E — — IC 


Let us calculate the contribution of the term E — — /C in eq. (jESj). The corre¬ 

sponding matrix element is: 


{cp\E - - t\dp') 

, fd f (K + 3/2)(K + 5/2)]' 

^ W- ? -i, 


Hp - P') 




(E7) 
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where Et is the intrinsic energy of the state t. Hence one obtains (cf eq. (jESj)): 




■p\E - / ‘‘P' 






(E8) 


2. Contribution due to the interaction between two protons in the continuum 


To illustrate the calculation of the interaction terms, we shall consider the Gaussian 
interaction: = Ho exp[—/3^(rj — YjY]. In the coordinate system T (cf sect. IIV Gil , 

the residual interaction between the two protons in the continuum is: 


y(res)^j^ _ sfiEf _ y^^-ffi{pcos{o)lyjl^f 

where fix is the dimensionless reduced mass of two protons. Hence, one obtains: 
(cp\VP“\A - 1, A)\c!p') = {K, p\VP“\A - 1, A)\K', p')S,,A,..pJ,,.rSLMSs.S'h,,.Ji 


(E9) 




= J da cos^(a) sin^(a)T^’^^(a)T'^r“(a) 

6{p - p') 


I'X. 


(p cos(o)/ 


X 


P 


5 St,t'Sl^,l'Jly,l'jL,L'Ss,S'Sj2p,J^p = ^c7''\p) 




(ElO) 


The Goulomb interaction between two protons in the continuum V^^\p) is calculated simi¬ 
larly. Hence, the hnal expression is (cf eqs. ()E5fl . (lElOfl i: 


lf.n(cp\Vf’-^‘\A-l,A) + VPE'(A-l,A)\Y, / dp'p'i\c'p')ujy\p') 

d 

E [drG)+d?>(p)] yy(p) 


(Ell) 


3. Contributions from the one-body potential 

Let us calculate the matrix element of the one-body potential uq between channels c and 
d in eq. (|E5|): 


{cp|no(H - 1) + vo{A)\cp) 


(E12) 
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Vo is the finite depth potential of a Woods-Saxon type with the spin-orbit term and the 
Coulomb potential of the daughter nucleus. Channels are defined in the coordinate system 
T. This system is not convenient for the calculation of the terms jnii) because Vq is a 
function of variables xi and X 2 which are associated with the coordinate systems Y (cf sect. 
IIV Cl aud Fig. [T]). The change from the coordinate system T to a coordinate system Y is 
done using the Raynal-Revai transformation 0 

Let us consider the transformation from the coordinate system (3) to the coordinate 
system (2) (cf Fig. P). The hyperspherical functions are transformed as follows: 


= E (l.A,llMKLyKl'‘'Hal) (E13) 

1x2 ’^2 

where (L 2 ^y 2 i^x 3 ^y 3 }KL are the Raynal-Revai coefficients. The hypermoment K and the an¬ 
gular momentum L are conserved by this transformation which corresponds to a rotation in 
the space of Jacobi coordinates. The summation over angular momenta 1^2 and ly^ associ¬ 
ated with directions X 2 and y 2 in the system Y defined by the coordinates (2) (cf Fig. [T]), is 
constrained as follows: 


Vx2~ ^y2\ ^ L < 1x2 + ^2 (E14) 

K = 2n -|- 1 x 2 ^j /2 


where n > 0 is an integer number. Analytical expression for the Raynal-Revai coefficients 
can be found in j^. Hence, one obtains: 


{cp\vo{A - l)\cp) 

= {t, K, {lx, ly)L, S, J2p, Jp\vo{A - 1)|F, K\ {I'x, l'y)L', S', J',; Jp') 

~ 'y ^ {^X2-)^y2\^x-}ly)KL{lx2'’^y2\^xi^y) 

111']' 

1x2 ’^^2 ’^^2 

X {K, {1x2, ly 2 )L, s, J 2 y, Jp\vo{A - l)\K', (4, 4)L', 4; Jp') (E15) 

Performing the transformation from L — S io j — j coupling scheme, one obtains for the last 
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term in eq. ()E15f) : 


{K, ly,)L, 5, J 2 p; Jp\vo{A - l)\K\ (4, 4)L', 5', J',; Jp') 


LSjlJy^L'S'jlJy 

3^2 ’3y2 ’3y2 


y2 


1x2 ^J/2 

Si S2 S' 


X 


/' /' T' 

‘'X2 ^2 

Si S2 

i' i' v 

Jx2 Jy2 '^‘ip 


jx2 jy2 

> {K, {jx2dy2)^ J2p^p\Vo{A - l)\K', U'^^Jy^), J2p]p) 


(E16) 


where: 


{K, {jx2,jy2), J2p;p\vo{A - 1)\K', J2p]p) 


, - ,rMA-i)r ^^ip-P') 

^=- ^i=^2A2^3x2p'x^\2>i'y^^3y2Ay^ = ^cc' \P) -:- 




(E17) 


pO Z ' yz ''‘'^'•'yz pO 

In this expression p^ is the dimensionless reduced mass of a system [A-2]®p and 
{pcos{a)/y/JI^) is the radial part of the one-body potential Vq for angular momenta 
^X 2 ) jx 2 - Similarly, one proceeds to calculate the term for a second proton, and eq. (Eg 
takes a form: 


(cp\v„(A - 1) + vM)Wp) ^ {VS(P)} (ei8 ) 

Hence, the total contribution from the one-body potentials becomes (cf eqs. (EH)) (lElSfl i: 
pi(c,p| {'i;o(^ - 1) + no(dl)} I j dp'p'^dp')J^]'^{p') 

c' 

= E{K?w}dv°(p) 

d 

(E19) 
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4. Contributions due to the interaction between protons in the continuum and 


nucleons in the daughter nucleus 


Let us calculate now the contribution from the residual two-body interaction between 
protons in the continuum and nucleons in the nucleus A-2 (cf eq. ()E5|1 ): 

j>A-l 


pHc,p\ [ dpp^c'p')uj^'^J’^{p') 

r-' ^ 


(E20) 


i<j 

In the second quantization form, this term is: 

j>A-i „ 

pHc,p\ Y. dp’p'l\c'p'}wpp°(,p') 

i<j d d 


= p^{c,pi - Y 




aiaifla. 


a,p,j,5 } L“q“/3J L“7 


(a ^ disc < /3 ^ coni') 
(7 G disc < 5 G coni) 

r 


. 05 ]^} \Y [ 

n' ^ 


(E21) 


/ /| I / /\ (+),0/ /\ 

pAcp)^)J KP) 


(E22) 


where is the reduced matrix element coupled to L of the residual interaction. Here 

we take into account the projection operator T by constraining summation over s.p. states 
a, (3, 7, 6: 


a G disc < /3 G cont 
7 G disc < 6 E cont 


The creation and annihilation operators are coupled in such a way that the operators acting 
on (quasi-)bound states are separated from those acting on the continuum states. One 
obtains: 




J' 


jp r 
' h js r 
j' j' 0 


V. 



(E23) 


49 






Inserting this expression in eq. ()E22ll and applying the Wigner-Eckart theorem, one gets: 




T^r -^2 


(a G disc < /3 G cont) 
(■7 E disc < 5 G cont) 

r, j' 


Ja j/3 r 


If J2p J 

J 7 js r 

> < 

If J 2 P J ^ 

J' J' 0 


f f 0 


Aa-fV'Wt') 


X (i^, (/., ly)L, S, J2p, p| I (alds)-^' 11 / dpV'i |ir', (4, /;)L', 5', J',, 


(E24) 


Let us consider the operator: 


^q^/3,7,(5(®^®5)m' (E25) 

(q: G disc, /3 £ cont) 

(7 G disc, 6 ^ cont) 

r 

in (ESH). Inserting the completness relation twice, one obtains: 


(q: G disc, j3 G cont) 
(7 ^ disc, 6 ^ cont) 

r 


^i/9,7d 

(a e diac,ep > 0) j P ,j2,j S ,j2dl,J2,MlM2 

(7 e disc, 65 > 0 ) 


r 


dridr2dr[dr'^rlr[^rlr'^^\jp,ri, j2,r2, 


X 0 '/ 3 , r 1, j2, r2, , Ml I (a^05) M /1 , r;, j2, ^2, ^2, M2) (, ri, j2, ^^2 > da, M2 1 


(E26) 


where an index j stands for the quantum numbers I and j of a given proton , e.g. ja denotes 
quantum numbers h, j 2 - and are the single particle energies associated with the states 
fd and 5. To simplify the demonstration, let us consider the case of two particles in the 
continuum which do not correspond to the same ensemble of quantum numbers j, i.e. jp 
differs from ja and, similarly, js differs from Applying the Wigner-Eckart theorem in 
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()E26fl . one finds: 


(cx £ disc , e0 > 0) 
(7 G disc , e § > 0) 

r 


5 / dridr2dr\dr'^rlr\rlr'^\jp,ri,j2,r2,Ji,Mi) 

jl3,j2,j5,j2yh,J2,MlM2 


x{J2M2fM'\JiMi){jfj2j2 < 


J /3 J 2 -^1 

J7 J2 
J' 0 J' 


U0{ri)us{r[) 5{r2-r'^),. , , , 

^- {js,r^,32,r2,J2,M2\ 


r , r ' ■^^’■^2 ^2 

'1 '1 '2 


(E27) 


where and ugir'^) are the radial wave functions of states /3 and 5, respectively. 

We shall sum now over energies of the states jd and 8 in (jE27jl . The term which depends 
on these energies can be written as: 


/ dradndridr[rir[ua{ra)u0{n)u0{ri)us{r[) 


X 


(E28) 


where v^^--^{ra-, Th) are the angle integrated, unsymmetrized matrix elements of the residual 
interaction, and the phase 0 equals: 0 = + js — T. Using the completness relation for s.p. 

states: 

^ + 00 


5^«en,U,r.(^)We„,/,i,r,(r') + j dc (r)Me,Z,i,r. (O = 5(r - U) (E29) 

n 

one hnds for the hrst term of 


'0 


/ dradndridr[rir[ua{ra)vl0 -^^{ra,rb)u^{ra)[8{rb - ri) - y^^Ung{rb)ung{ri)] 

n/3 

x[ 8 {rb-r[) -'^Uns{rb)uns{r[)] (E30) 


ns 


where Unfairb) and Ung{rb) are the radial functions of proton (quasi-)bound states with quan¬ 
tum numbers = [// 3 , j/ 3 ] and js = [h,js], respectively. Similarly, the second term in (IE28fl 
becomes: 


-1)'^ / dradrbdridr[rir[ua{ra)vl 0 ^^-{ra,rb)uj{rb)[ 8 {rb - n) - M„^(rb)M„^(ri)] 

np 

x[(5(r„ - r[) - ^Ung{ra)ung{r'^)] (E31) 


ns 
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Using expressions ()E27fl . ()E28fl . (jESO)), and (jEdlfl . the operator ()E25jl can be written as: 




(a £ disc, j3 ^ cont) 
(7 G disc, 5 G cont) 

r 


E E f{J'fj2W2M2J'M'\JiMi) 

(a e disc) j/3d2,js,Jl,J2,MlM2 


(7 G disc) 

r 


J/3 J2 

J7 72 J2 
J' 0 J' 


dradndridr[rir[dr 2 rl x \jf3,ri,j2,r2, Ji, Mi){js,r[,j 2 ,r 2 , J 2 ,M 2 \ua{ra 


X i j(r„,rb)[(5(rb - n) - Mn^(r6)M„^(ri)][(5(rb - UJ -J2^^sin)uns{r[)] 


np 


ns 


-{-l)\{rh)vl^^ -s -{ra,rb)[6{rh - ri) - Mn,(ra)Mn,(r’i)] 


n /3 


715 


(E32) 


or as: 


E E f (J')""-^ 2/2 (J 2 M 2 J'M ViMi) 

(a e disc) jpd2,js,Jl,J2,MlM2 
(7 G disc) 

r 


jp 

^ h 
J' 


32 Jl 

32 J 2 
0 J' 


dr2r2dridr[rir[ 


I j/3, r-i, j 2 , ^ 2 , Jl, Ml) /(r 1 , r'l) {js , r'l, js, r 2 , J 2 , M 2 


(E33) 


The form of the operator firi,r[) in the above expression can be found easily from ()E32|1 . 
Let us now calculate the operator (TeH^ for antisymmetrized two-particle states \ah) and 
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led): 


(-^1 E E 

(cediso) j/3j2,jsJl,J2,MlM2 


(7 E disc) 

r 


J/3 J2 Jl 

j'y 32 J 2 
J' 0 J' 


> / dr2rldridr[rir[ 


X \j/3, rij2, r2, Jl, Mi)/(ri, r^, js, ?^2, J 2 , Malcd) 

= E E T{jfJ232{J2M2J'M'\J,Mi) 

(a e disc) jl3,j2,js,Jl,J2,MlM2 
(7 E disc) 

r 

x(a6|j/3,ri,j2,r2, Ji,Mi)(„a 5 ) 2f{ri,r[) {nas){js,r[j 2 ,r 2 , J 2 ,M 2 \cd) 


jp 32 Jl 
J 7 d 2 J2 
J' 0 J' 


> / dr2r2dridr[rir[ 


(E34) 


where 'nas' means that the considered state is not antisymmetrized. Formally, above 
expression suggests that the operator (lE^ acts only on the hrst proton in different 
unsymmetrized states referenced by ri and r[. Indeed, the particle in a state |j 2 ,''" 2 ) enters 
in the calculation only through the angular momentum coupling. 


To calculate the last term in ()E24fl . one goes from L-S to j-j coupling scheme and 
then applies the Raynal-Revai transformation to obtain: 


(A', (A. l,)L, S, .h„ p| |(ay-j)"'ll A', (4,4)L', S', 4 , p') 

= E E LSLivd'S'iijl 

1x2 3y2 3x2 3y2 3x2 dy2 ’3x2 '3y2 


X 


^X2 

ly2 

L 


' ^'x2 

V 

L' 

Si 

S2 

S 

> < 

Si 

S2 

S' 

3x2 

3y2 

J2p 

} 


3x2 

dy2 

J' 

'J2p 


^ {K, {3X2,3y2), J2p, p\\{ald;s)'^'\\K', ^p, p) 


(E35) 


As discussed above (cf expression (IE34fl h one may consider that only a proton labeled by 
jx 2 on the l.h.s., and a proton labeled by on the r.h.s. enter in this calculation. Hence, 
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the reduced matrix element becomes: 


(■^5 Ux2i 3y2)i P\ I iflfjOjs) I \K , (jx2j Jj/ 2 )’ '^2p5 P ) 

jx2 jy2 '^‘2p 
i' i' V 

Jx2 Jy2 2p 

J' 0 J' 


J2pJ 2p{J ) jy2 ^ 


> I dada cos^(a) sm^(Q;) cos^(q;') sin^(Q;'' 


^I^2 3y2( m^-2’^y2( >) Up{p COS a/^) usjp' cosa' / ^) 6{psma - p' sing') 
^ p cos a/y/JI^ p'cosa'/y^ (psina)^ 

X ^ 1.2 h ^j.2 Pp ^i'^2 ^^-2 P^ ^^y2 Py2 


(E36) 


Using this last expression and (jEddjl one can write (lE^ as: 


-pV2 




'■ + CXD 


|^(r, a, p, 7, <5, ^2, p) + B{T, a, p, 7, 6, ^2, p)} 


c', (a, 'y) ^ disc, E, j' 

^y2 ’ ^y2 ’ ^P ’ ^0 ’ ^<5 > 


.2 , - ' 7 I do j'y J [ J It It H 


xr^jd(_i).'^ " W - - " ^(t||(atd,K||t') 

36 dp r 1 I j'p hp J 


^ Px2py2VxPy) KLpx2i^y2\^xi^y) K' L' LSj pj^L S jy^ J2pJ2p 


X < 


Ip ly2 L 

^ Si S2 S 

> < 

< ' 

^5 ^y2 ^ 

Si S 2 S' 

ij 

^ J2p J'2p J' \ 

dp dp2 ■1‘ip 

V / 


js jy2 J 2 P 


^ 3s dp dj/2 j 


where 0' = 1 + + jy^ 


i2p -r ^2p 


r + /;. 


(E37) 


A(r, a, P, 7, 6, ly 2 , p) in (IE37|) is equal to: 

A{T,a,P,'-f, 5 ,ly^,p) = j dradrfeMa(ra)Mp(rb)nr^-j(ra,rb)M^(ra)M 5 (rb) 

X J dp'p'idada cos^ a sin^ a cos^ a sin^ (q;)\ 1/^’/^^ (a') 

up{p cos a/y/pP)us{p' cos a'/y/JIp) 5{psm. a - p' sm. a') (+), 0 / 


X- 


p cos a Z^/pP p'cos a'/ y/p^ 


(p sin a) 2 




(E38) 
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and -B(r, a, /3, 7 , S, p) is: 


B(T,a,(3,-f,6,ly2,p) = J dradrbUa{ra)up{rb)vl0^{ra,rb)u^{rb)u5{ra) 

X J dp pi dada cos^ a sin^ a cos^ a sin^ (a') 

up{p cos a/y/JI^) us{p'cos a'/y/p^) 6{p sin a - p'sin a') (+), 0 / , 


X 


pcosa/y/p^ p'cosa'/y/p^ 


(psin a )2 




(E39) 


Using the completness relation for s.p. states and integrating over energies e /3 and es, the 
hrst term in ()F;37j) becomes: 


^+00 


/ de,de,A{T,a,f5,-f,5,ly^,p) 


= / dp'pU dvadvbdada' cos^ a sin^ a cos^ a' sin^ a 


2 „/„• 2 „/ 




XMa(ra)ni^^j(ra, (a') 


X 


(5(psin a — p' sin a') 
(psino;)^ 


'5{rb- pcosa/^/p^) Un^ {rb)un 0 (p cos a/^/JI^) 


n 


{pcosa/y/jQ 


X 


5{rb - p'cos a'/y/JI^) Uns{rb)uns{p' cos a'/y/JI^) 


n 


(p' cos a'/,/p^) 


(+), 0 / /\ 

sV (d) 


(E40) 


The role of terms: 


Un 0 {rb)un 0 (p COS a/^/p^)/{p cos a/y/pZ) 


and 

Uns {n)Uns (P' COS tt'/ y/p^) / {p COS ol jy/JQ 


in (irail is to project on the (qnasi-)bonnd s.p. states. Their ontcome will be taken into 
acconnt by the method described in appendix |Fl and, at this point, we remove them. Inte¬ 
grating over rb, p' and 0 : 2 , one can rewrite (jEdOfl as follows: 


da cos^ a sin^ (a) / draUo,{ra)u^{ra)v'^- 0 ^b{ra, p cos a/ (a)a;^^^/’''(p) 


Is ,ly 




(E41) 


Similarly, nsing the completness relation for s.p. states and integrating over energies ep 
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and e^, the second term in (IE37fl takes the form: 


^+oo 


P ^ I o;, /3, 5, ly2-, p) 


— _ ( _ 


I dp p dradr^dada' cos^ a sin^ a cos^ a' sin"* a 


Jl „/ „• 2 „/ 




^(nsin a — p' sin a') 

^ - r ~-—^ 

(psm ay 


6{rb- p cos a/^/JQ (n)Mn^ (p cos a/yjQ 


n 


X 


5{ra- p' cos a'/y/jQ Uns{ra)Uns{,p' cosa'/y/jQ 


{pcosa/y%) 


(E42) 


Ta {p' cosa'/yjQ 

Integrating over r^, and a, and neglecting the projection operators on the (qnasi-)bonnd 
s.p. states, one obtains: 

— (—j dp p^ p~^ da' cos a' sin^ a' 

XM„(p' cosa /^/py)v]^^-^-{p' cos a'/ pcosa°/^/p;^)n^(pcosQ;°/ 

where: 


n . / p sma2 

a = arcsm 

P 


Using expressions (jEill) and (lE^ . one can calculate the term (lE^ . Let us write this 
term in the form: 


E [hhU)+U”'b)] ‘^y\p) 


(E44) 


where V}lf''\p) and v}^’‘\p) are, respectively, local and non-local potentials. Using expres¬ 
sions (lE8l) . (jllH), (jHni), jEHlj) and (IE44I1 . one can write eq. dlED as: 


E-E, 


(i^ + 3/2)(iL + 5/2) j' 


2 m 1^ df)- 


P" 




or: 


E (U”‘’(p) + hfU) + k“’(p) + U“'(p) + vkU)) dt’’°(p)= 'pM 


(E45) 


E 


Hcc'{p)^^p]'^{p) = Wj^yp) 


(E46) 
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where the channel-channel conpling potentials are: 


Hcc\p) — 


E-E, 


+K 


r_^ 

2m \ dp^ 

+ v^:;\p)+v^:^{p) + v;:i?^\p)+v^\p) 


{K + ?,/2){K + b/2) \ 


(E47) 


APPENDIX F: PROJECTION OPERATOR IN T-SUBSPACE 

In appendix lEl we have discnssed the calcnlation of the fnnction neglecting the 

projection operator T: 

= (FI) 

Below, we shall show how this operator can be inclnded effectively, i.e how one can solve 
eqnation: {E — = \wj). By dehnition, belongs to T. Let ns now dehne 

’forbidden states’ for the state In the coordinate system Y (cf sect. II V Cfl . they can 

be expressed as: 


^t,n,Jt,n,s,J fjn) y y js^ J) (^2) 

Here, t{^) is a state of the nnclens A-2 and Un{xi) is a (qnasi-)bonnd state of the one- 
body potential Hq. It is the total angnlar momentnm of a state t{^) and is the total 
angnlar momentnm of the state u^^ixi). It and jn are conpled to Jt,n- {fsiVi)} in (jF2jl is a 
complete set of spline fnnctions and js is the angnlar momentnm associated with the 
spline fnnction fg. Jt,n and js are conpled to the total angnlar momentnm J. To simplify 
notation, we shall write \^t,n,Jt,n,s,j) as |‘he,j), where e = {t,n, Jt^nys). States |$e,j) satisfy: 

T|d>e,j) = 0 (F3) 

and are called “forbidden” in the sense that as belongs to T then: 

= 0 (F4) 

for all e. The set of states |‘he,j) is constrncted from all considered states |t) of a nnclens 
A-2 and all (qnasi)bonnd states of a potential hg, nsing the complete set of spline fnnctions 

{fsiVi)}- 
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Solution of eq. (gni) is written in the form: 


|u,W> = |^W'“) + ^A.|u.(+)) (F5) 

e 

i.e. we want to determine coefficients Ag. ) in eq. is a solution of the equation: 

(B-/f)l4+>> = |tg,o> (F6) 

The calculation of |aJe^*) is analogous to the calculation of (cf appendix El- From 

CD and (iFl one can see that fulhlls the equation: 

(B-//)|w]+>) = |u,,) + ^A.|4g,r) (F7) 

e 

Coefficients Ae in (jF5|) are chosen to assure that belongs to T, i.e. they satisfy: 

{a;f)|$e,j) = 0 Ve (F8) 

or equivalently: 

=0 V 6 (F9) 

e' 

Hence, belongs to T and is the solution of equation (cf eq. IF7fl : 

(B-//f)|wy>) = l«’r) + E^«l«'»> (FlO) 

e 

Multiplying eq. (IFlOf) from the l.h.s. by the operator T, one obtains eq. (14(1 II . what completes 
the argument. 
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